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ABSTRACT 

T“ 

The  1988  STLE  Annual  Meeting  held  a  two- 
session  symposium  on  cavitation  in  fluid  films. 
This  paper  Is  an  edited  review  of  the  current 
research  as  discussed  by  the  invited  participants. 
Phenomena  and  experimental  observations  Included 
gaseous  cavitation,  vapor  cavitation,  and  gas 
entrainment.  Cavitation  in  flooded,  starved,  and 
dynamically  loaded  journal  bearings,  as  well  as 
squeeze-film  dampers  was  reviewed.  Observations 
of  cavitation  damage  In  bearings  and  the  possi¬ 
bility  of  cavitation  between  parallel  plates  with 
microasperities  were  discussed.  The  session  on 
theoretical  modeling  summarized  the  transcavity 
fluid  transport  process,  meniscus  motion,  and 
geometry,  or  form,  of  the  film  during  rupture  and 
reformation.  Performance  effects  were  related  to 
heat-transfer  models  In  the  cavitated  region  and 
hysteresis  Influences  on  rotor-dynamic  coeffi¬ 
cients.  .Contributors  presented  a  number  of 
cavltatlotwlgorlthms  together  with  solution 
procedures  using  the  finite-difference  and  finite- 
element  methods.  Although  Newtonian  fluids  were 
assumed  in  most  of  the  discussions,  the  effect  of 
non-Newtonian  Ifluids  on  cavitation  was  also 
discussed.  \ 

INTRODUCTION  eft)  ^ - 


This  paper  Is  a  survey  of  the  state-of-the- 
art  In  research  in  cavitating  fluid-film  bearings 
and  seals  based  on  a  symposium  held  at  the  STLE 
Annual  Conference  In  Cleveland,  Ohio,  in  May 
1988.  The  purpose  of  this  article  is  to  discuss 
the  fundamental  Issues,  experimental  observations, 
theoretical  modeling  techniques,  and  further 
research  needed  In  this  field.  It  is  hoped  that 
this  paper  can  serve  as  a  reference  publication 
for  both  the  experienced  researcher  and  those  who 
are  becoming  Interested  in  this  area.  Neverthe¬ 
less,  this  article  Is  not  a  general  review  paper, 
as  it  is  largely  based  on  the  viewpoints  of  the 
participants  In  the  symposium.  Each  participant 
provided  a  short  position  paper  covering  his 
presentation  from  which  the  editors  have  drawn 
freely. 


This  article  is  divided  into  three  parts.  In 
the  first  part,  we  present  experimental  results  to 
provide  the  reader  with  the  observed  phenomena  as 
well  as  some  of  the  pertinent  measurements;  the 
second  part  of  the  paper  deals  with  theoretical 
modeling  techniques.  Part  III  consists  of  the 
editor's  summary. 

In  the  text  that  follows,  different  types  of 
cavitation  are  described: 

Gaseous  Cavitation 


The  lubricating  oil  often  contains  some  dis¬ 
solved  air  and/or  other  gases.  It  is,  therefore, 
likely  that  these  gases  could  come  out  of  the 
solution  If  the  fluid  pressure  falls  below  the 
saturation  pressure  of  the  dissolved  gases  (nearly 
atmospheric  pressure).  This  type  of  cavitation 
is  called  gaseous  cavitation.  In  journal  bearings 
under  steady-state  conditions,  one  would  generally 
expect  gaseous  cavitation  to  take  place  In  the 
divergent  region  where  the  mating  surfaces  tend  to 
pull  the  lubricant  apart,  thus  creating  the  condi¬ 
tion  in  which  the  gases  could  emerge  from  the 
solution  and  give  rise  to  "atmospheric  pressure". 
Side  entrainment  is  another  manifestation  of 
gaseous  cavitation  when  bubbles  at  the  ambient 
pressure  of  the  sides  are  drawn  into  the  negative 
pressure  regions  of  the  bearing. 

Vapor  Cavitation 

If  the  local  pressure  of  the  fluid  coincides 
with  that  of  the  vapor  pressure,  then  the  boiling 
of  the  fluid  by  pressure  reduction  causes  the 
formation  of  a  so-called  vapor  cavitation.  This 
type  of  cavitation  normally  occurs  in  dynamically 
loaded  bearings  used  in  common  machinery. 

Understanding  the  mechanism  of  the  formation 
and  collapse  of  cavitation  bubbles  Is  of  funda¬ 
mental  Importance  since  it  Is  known  to  affect  the 
performance  of  the  bearing.  Furthermore,  it  is 
the  vapor  cavitation  that  is  known  to  be  harmful 
as  the  bubbles  violently  collapse  against  the 
bushing  surfaces,  resulting  In  cavitation  damage. 
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PART  III  —  EDITORS  SUMMARY 
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PART  I 

FUNDAMENTAL  AND  EXPERIMENTAL  OBSERVATION 
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NOMENCLATURE 


b  Pad  width,  m 

C  8earing  radial  clearance,  m 

D  Bearing  diameter,  m 

d  Depth  of  ol 1 ,  m 

e  Eccentricity,  m 

F  Force,  N 

f  Frequency  of  load  (Hz),  s-' 

G  Duty  parameter,  CpUb/H] 

h  Film  thickness,  m 

ha  Amplitude  of  oscillation,  m 

L  Length  of  bearing  In  z-directlon,  m 
N  Shaft  rotational  frequency,  rev/mln 

P  Unit  loading  (N/LD) 

ps  Oil  supply  pressure,  m3/s 

Q  Volumetric  oil  flow,  m3/s 

Q  Q/(0.5  UCL) 

R  Bearing  or  journal  radius,  m 
t  Time,  s 

U  Linear  velocity,  m/s 

u  Velocity  of  the  fluid  In  the  circumfer¬ 
ential  direction,  m/s 

N  Load  on  bearing,  N 

x  Coordinate  In  circumferential  direction,  m 

y  Coordinate  in  radial  or  normal  direction,  m 


y0  Initial  film  thickness,  m 

z  Coordinate  in  axial  direction,  m 

(i  Fluid  film  arc,  rad 

e  Eccentricity  ratio  of  bearing  (e/c) 

0  Coordinate  in  plane  of  rotation,  rad 
60  Angle  of  line  of  centers,  rad 

P  Viscosity,  Ns/m2 

Pt  Coefficient  of  friction 

a  Shear  stress,  N/m2 

t  Period  of  dynamic  load,  s 

4>  Attitude  angle,  (0O  -  «•),  rad 

w  Angular  velocity  of  shaft,  rad/s 

Subscripts: 

B  Beginning  of  bearing  pad 

d  Dynamic 

E  End  of  bearing  pad 

F  Ful 1  fluid  fl 1m 

min  Minimum 

S  Starvation 

s  Static,  supply 

st  Steady 

1  Start  of  fluid  film 

2  End  of  fluid  film 
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EXPERIMENTAL  STUDIES  IN  CAVITATION 
H.  HESHMAT' 

The  experiments  cited  here  deal  primarily 
with  the  flow  and  formation  of  oil  streamlets  In 
cavltated  regions. 

FLOODED  BEARINGS 


Consider  Fig.  1,  which  depicts  a  hydro- 
dynamic  journal  bearing.  As  shown,  the  fluid 
film  ends  at  ©2  after  which  it  cavitates; 
the  space  at  0  >  ©2  is  filled  with  oil  stream¬ 
lets,  released  gasses,  ambient  air,  and  possibly 
oil  vapor.  In  the  absence  of  hydrodynamic  posi¬ 
tive  pressures  In  the  region  ©2  -  ©£.  the  amount 
of  flow  (©2)  leaving  Q2  remains  constant  at 
all  ©  >  ©2-  If  one  observes  the  cavitation 
region  through  the  bearing  t  ansparent  shell,  one 
sees  fluid  film  streamlets  in  the  manner  shown  In 
view  B  of  Fig.  1.  A  sample  of  such  experimen¬ 
tally  observed  streamlets  seen  through  a  transpa¬ 
rent  bearing  shell  Is  shown  In  Fig.  2  (cf.  (1>). 
This  figure  reveals  that  the  oil  adhering  to  the 
stationary  surface  is  in  narrow  strips.  The 
fluid  film  emerging  from  a  cavltated  region 
beyond  the  physical  end  of  the  pad  ©£  is  shown 
in  Fig.  3.  The  pattern  of  ridges  and  troughs  in¬ 
dicates  that  next  to  the  shaft  is  a  sublayer  of 
fluid  that  completely  covers  the  rotating  shaft. 
This  is  unlike  the  stationary  surface.  The  fluid 
as  it  emerges  from  the  bearing  pad  does  not  leave 
the  shaft  but  adheres  to  it  even  when  the  film  is 
not  restrained  by  an  upper  bearing  surface  but  is 
open  to  the  atmosphere.  The  streamlets  thus  will 
adhere  to  the  shaft  even  when  they  cross  an  oil 
groove.  The  streamlets  emerge  at  the  top  of  the 
picture  from  a  cavltated  zone  and  cross  an  axial 
groove  before  entering  the  downstream  pad.  As 
seen,  neither  the  axial  distribution  nor  the  indi¬ 
vidual  width  of  the  streamlets  changes  in  cros¬ 
sing  the  groove,  thus  indicating  adhesion  of  the 
entire  flow  (Q2>  to  the  moving  surface. 

The  above  discussion  suggests  that  the  shape 
of  the  fluid  film  in  the  cavitated  region  is  in 
the  form  depicted  in  Fig.  4.  The  fluid  adheres 
In  thin  strips  to  the  stationary  surface,  but 
forms  a  continuous  sublayer  on  the  moving  surface. 

This  pattern  prevails  as  long  as  there  is  an 
upper  pad  and  a  thin  spacing  gap  between  shaft 
and  bearing  surface.  Upon  termination  of  the  pad, 
the  tips  detach,  and  the  entire  flow  (Q2)  adheres 
to  the  shaft  in  the  manner  shown  in  Fig.  3.  This 
"detached  pattern"  also  prevails  in  the  oil 
groove,  whose  dimensions  are  two  to  three  orders 
of  magnitude  larger  than  the  thickness  of  the  film 
layer  and  is  thus  equivalent  to  an  open  space. 

The  detached  pattern  reveals  that  there  is  a  very 
weak  bond  between  moving  streamlets  and 
stationary  surface,  permitting  side  entrainment 
of  gas  to  take  place  in  cavitated  regions. 
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Fig.  1.  -  Cavitation  in  journal  bearings. 


TEMPERATURE  EFFECTS 


The  next  topic  pertains  to  the  temperatures 
of  the  hot  layer  (Q2).  Numerous  tests  have  shown 
that  the  temperatures  of  the  fluid  film  in  the 
cavitated  region  are  more  or  less  constant.  This 
is  true  despite  seemingly  further  energy  dissipa¬ 
tion,  which  takes  place  due  to  viscous  shear  over 
the  extent  of  the  cavltated  region.  Figure  5 
shows,  for  five  sets  of  experimental  operating 
conditions,  that  the  temperatures  of  the  cavl¬ 
tated  region  remained  constant  even  though  this 
region  extended  over  nearly  80®.  A  possible 
explanation  for  this  phenomenon  is  that  the 
entrainment  of  ambient  air  separates  the  fluid 
from  the  bearing  surface  and  causes  less  heat 
conduction  to  the  thermocouples;  i.e.,  thermo¬ 
couple  tips  are  exposed  to  gas  rather  than  hot 
oil.  A  partial  explanation  has  also  been  given 
in  terms  of  released  gases  and  oil  vaporization, 
which  requires  high  enthalpies.  However,  even  in 
cases  where  the  subambient  pressure  is  not  low 
enough  to  cause  oil  vaporization,  a  constant 
temperature  has  been  observed.  To  this  end,  the 
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L/D  =  1;  N  =  1000  rpm;  P  =  60  psi 
(C/R)  =  2x  10-3  ;  Ps  =  0.5  psi 

2.  -  Flow  of  oil  streamlets  over  groove. 


View  A-A 

Fig.  4.  -  Oil  streamlets  In  cavltated  regions  or 
journal  tarings. 
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Fig.  5.  -  Constant  temperatures  In  cavltated  re 
glons;  P  =  2.76  MPa. 


shape  of  the  fluid  film  as  postulated  in  Fig.  4 
offers  an  explanation  of  the  role  of  gas  entrain¬ 
ment  and  constancy  of  temperature  in  the  cavitated 
region.  The  lower  part  of  Fig.  4  shows  the  trans¬ 
verse  velocity  distribution  in  the  oil  streamlets. 
As  seen,  most  of  the  oil  moves  with  a  velocity 
very  near  the  journal  speed  <U);  the  principal 
strain  thus  occurs  near  the  stationary  surface. 
Since  the  power  dissipation  is  related  to  the 
stress  o  =»  o(du/dy)  at  the  moving  surface  and 
since,  here  (du/dy)  =  0,  the  losses  would  be  neg¬ 
ligible.  Consequently,  there  would  be  no  rise  in 
temperature. 

To  show  that  oil  generally  adheres  to  the 
moving  surface,  a  series  of  tests  was  conducted 
with  a  tapered-land  thrust  bearing  with  reverse 
rotation.  The  apparatus  and  the  experimental  pro¬ 
cedures  are  described  by  Heshmat  <2).  The  test 
bearing  was  instrumented  with  pressure  taps  and 
thermocouples  to  monitor  the  fluid  film  pressure 
as  well  as  temperature.  Rotation  was  opposite  to 
the  customary  direction  so  that  only  the  flat  por¬ 
tion  was  loaded.  Figure  6  gives  some  of  the 
results  obtained  with  this  series  of  tests.  The 
noteworthy  points  pertinent  to  the  cavitation  ate 
the  following: 

1.  The  pressure  cap  (located  at  the  center 
of  the  flat  portion  of  the  pad)  response  attests 
to  the  presence  of  a  hydrodynamic  film,  not 
boundary  lubrication  over  the  pad  at  the  onset. 

2.  In  the  diverging  portion  of  the  film,  the 
temperatures  actually  tend  to  drop  (Fig.  6). 

This  further  supports  the  view  that,  in  the 
diverging  film  with  cavitation,  temperature  rise 
is  negligible  because  of  the  negligible  power 
losses  mentioned  earlier.  Furthermore,  the  end 
flow  and  streamlets  in  the  cavitation  region  - 
induced  by  the  diverging  wedge  configuration  - 
basically  adhere  to  the  moving  surface  indepen¬ 
dent  of  surface  geometry  (journal  or  flat  collar). 


Fig.  6.  -  Temperature  profiles  in  a  flat,  diver¬ 
gent  bearing. 


STARVED  BEARINGS 

An  experimental  study  documenting  the  perfor¬ 
mance  of  journal  bearings  for  conditions  ranging 
from  full  hydrodynamic  lubrication  to  states  of 
extreme  lubricant  starvation  is  described  by 
Heshmat  and  Pinkus  (3).  The  effects  of  an  incom¬ 
plete  fluid  film  on  the  extent  of  the  film  and 
the  formation  of  cavitation  will  be  reviewed  here 
as  functions  of  the  amount  of  lubricant  available 
to  the  bearing. 

Figure  7  presents  the  hydrodynamics  of  a 
starved  bearing  vis-a-vis  a  full  film  or  flooded 
condition.  For  a  fixed  speed  and  supply  oil 
temperature,  the  deviations  of  a  starved  bearing 
from  one  operating  with  a  full  fluid  film  are  as 
fol lows: 

1.  The  film  starts  later  and  terminates 
earlier,  i.e,  0i  >  0]p  and  ©2  <  ©2F>  producing, 
in  essence,  something  similar  to  a  partial 
bearing  -  though  the  upstream,  and  sometimes  also 
the  downstream,  boundary  conditions  are  different. 

2.  The  decrease  in  the  fluid  film  extent  in 
starved  bearings  is  shown  in  Fig.  8  and  is  pretty 
close  to  the  theoretical  values  reported  else¬ 
where.  The  start  and  end  of  the  film  was  ascer¬ 
tained  by  means  of  pressure  probes  which  located 
the  points  of  the  commencement  and  termination  of 
the  hydrodynamic  pressures. 

3.  The  behavior  of  the  fluid  film  in  the 
starvation  zone  (0g  -  0j)  was  observed  via  a 
transparent  bearing  shell.  Its  appearance  is 
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Fig.  7.  -  Full  and  starved  fluid  films. 
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a.  Basic  Film  Pattern 


/  N  =  1800  rpm 

/  SAE  20  Oil 
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Q2 

Fig.  8.  -  Extent  of  fluid  film  In  starved 
bearings. 


sketched  In  Fig.  9.  As  seen,  the  full  film  under 
starved  conditions  started  at  more  or  less  a 
constant  angle,  0i  »  ©i p  and  did  not  vary  along 
z.  The  starved  zone  between  the  start  of  the  pad 
and  ©i p  underwent  cyclic  changes  in 
appearance.  Most  of  the  time,  the  zone  0g  S  9  S 
01  was  filled  with  an  amorphous  mixture  of  oil, 
foam,  and  gas  bubbles  (Fig.  9(a)).  However, 
every  0.5  to  1.0  s,  a  pulse,  originating  from  the 
cavitated  zone  downstream  (0  >  0%) ,  transformed 
this  region  into  a  pattern  of  oil  streamlets 
(Fig.  9(b))  similar  to  the  conventional 
downstream  cavitation.  The  pulse  frequency  was 
checked  against  the  journal's  rotational 
frequency,  but  no  relation  nor  connection  was 
found  between  them.  However,  the  pulsating 
streamlets  and  their  frequency  were  related  to 
the  degree  of  starvation.  The  higher  the  degree 
of  starvation,  the  longer  the  period  of  pulse 
frequency. 

4.  The  upstream  portion  of  the  starved 
bearing  0g  §  0  i  ©i  filled  with  an  oil  and  gas 
mixture  is  shown  schematically  in  Fig.  9(a). 
Although  gaseous  bubbles  were  visible  through  the 
transparent  bearing,  they  disappeared  in  the  film 
and  had  no  effect  on  the  streamlets'  patterns  in 
the  cavitated  zone. 

Based  on  these  observations,  the  model 
provides  a  reasonably  accurate  prediction  as 
given  in  the  literature  (1,4-6). 


b.  Pattern  During  Periodic  Pulses 
Fig.  9.  -  Fluid  film  pattern  in  starved  bearings. 


DYNAMICALLY  LOADED  BEARINGS 

Film  rupture  in  dynamically  loaded  bearings 
has  been  experimentally  observed  (7-24)  to  occur 
in  two  general  patterns:  (1)  film  rupture  always 
exists  through  the  load  cycle,  and  (2)  occurrence 
and  disappearance  of  the  film  rupture  are  repeated 
every  load  cycle.  In  the  case  of  a  nutating 
squeeze-film  damper,  the  ruptured  film  travels 
circumferentially  ahead  of  the  minimum  film 
thickness  (25). 

Recent  studies  by  Nakahara  et  al .  (18), 

Nakal  et  al.  (14),  and  Natsumeda  and  Someya  (26), 
all  of  whom  used  special  pressure  transducers 
(semiconductor  strain  gauge  type)  mounted  in  the 
test  journal  aid  bearing,  report  the  existence  of 
a  sharp  negative  pressure  except  under  a  lightly 
loaded  condition.  Negative  pressure  (tensile)  up 
to  -1.2  MPa,  as  shown  in  Fig.  10,  was  reported  by 
Natsumeda  and  Someya  et  al .  (26).  It  is  concluded 
that  the  heavier  the  dynamic  load,  the  greater  the 
negative  pressure.  Such  a  large  value  has  been 
scarcely  reported  by  others.  Dyer  and  Reason 
(11)  reported  tensile  stress  with  a  maximum  value 
of  0,64  MPa.  However,  their  experimental  results 
showed  that  the  magnitude  of  tensile  stress 
decreases  as  the  eccentricity  ratio  is  increased. 
Nakahara  et  al .  (18)  concluded  from  their  experi¬ 
ments  that  cavity  pattern  changes  were  caused  by 
the  cavity  growth  due  to  the  diffusion  of  dis¬ 
solved  gas  in  oil  under  constant  speed.  Further¬ 
more,  as  reported  recently  by  Walton  et  al .  (25), 


Fig.  10.  -  Circumferential  pressure  distribution 
of  journal  bearing  under  dynamic  load  (26). 


cavity  patterns  changed  with  the  change  of  cavity 
locations.  Both  subcavity  pressure  and  radial 
flow  resistance  across  the  circumferential  groove 
changed  discontlnuously  (traveling  pressure  wave 
up  to  five  times  that  of  supplied  pressure). 


Depending  on  the  apparatus  and  test  operating 
conditions,  five  or  more  modes  of  transitional 
phenomena  with  cavity  change  ware  categorized  by 
the  Investigators.  Various  boundary  conditions 
have  been  used  to  predict  dynamically  loaded  bear¬ 
ing  performance.  These  are  Sommerfeld  conditions 
for  light  load,  Reynolds  conditions  for  heavy 
load,  separation  boundary  condition  for  detecting 
tensile  stress  in  the  oil-film  region  downstream 
of  the  minimum  film  thickness  position,  and 
others  In  order  to  suit  the  fulfillment  of  the 
experimental  data  pattern.  While  the  agreement 
between  analysis  and  test  is  of  the  right  order 
of  magnitude  (depending  on  the  test  and  boundary 
conditions  assumed)  In  predicting  the  extent  of 
the  film  rupture  zone,  analysis  does  not  ade¬ 
quately  reflect  the  noted  variations  in  the 
boundary  conditions  observed  and  does  not  account 
for  the  change  in  the  shape  of  the  film  rupture 
zone  with  Increasing  speed. 


In  a  review  of  experimental  results  described 
by  Walton  et  al.  (25),  pressure  variations  In  the 


damper  grooves  exhibited  significant  flow  rever¬ 
sals  (circumferential  feed  and  drain  grooves  in 
the  sealed  squeeze  film  damper).  Also,  the 
grooves,  instead  of  being  full  of  oil,  appeared 
to  experience  film  rupture,  which  contradicts  the 
boundary  condition  assumption  that  the  groove  acts 
as  a  large  sink.  The  flow  fluctuations  in  the 
damper  oil  grooves,  while  not  fully  explainable, 
do  affect  damper  performance  and  various  stages  of 
cavitation  formation. 

It  Is  believed  that  side  entrainment  of  gas 
takes  place  In  most  of  those  tribological  test 
apparatus  and  sets  up  various  cavitation  config¬ 
urations.  These  discrepancies  could  be  attri¬ 
buted,  at  least  partially,  to  the  differences  In 
the  interpretation  of  the  various  observations. 


CONCLUDING  REMARKS 


Based  on  the  review  of  the  literature  and 
results  of  an  experimental  Investigation,  the 
following  remarks  are  made: 

1.  Side  entrainment  of  the  surrounding  gas 
to  the  bearing  gap  at  the  downstream  position  Is 
conceivable. 

2.  The  shape  of  the  fluid  film  in  the  cavi- 
tated  region  Is  in  the  form  given  by  Fig.  4.  The 
fluid  adheres  by  means  of  weak  thin  strips  to  the 
stationary  surface,  but  forms  a  continuous  sub¬ 
layer  on  the  moving  surface.  This  pattern  pre¬ 
vails  as  long  as  there  is  an  upper  pad.  Upon 
termination  of  the  pad,  the  tips  detach,  and  the 
entire  exit  flow  adheres  to  the  moving  surface. 

In  this  region  there  would  be  no  rise  In 
temperature. 

3.  Negative  pressure  (tensile  stress)  has 
been  detected  experimentally  on  the  order  of 

1 .2  MPa  (absolute). 

4.  To  fully  understand  and  explain  a  sharp 
dent  of  negative  pressure,  one  needs  to  consider 
the  lubricant  as  a  two-phase  liquid  consisting  of 
oil  and  bubbles  of  gas  and  to  take  into  account 
the  surface-dilation  effect  of  the  surface  ten¬ 
sion  together  with  tensile  strength  of  the  film 
in  the  governing  equation  of  motion. 

5.  Under  dynamic  conditions,  other  types  of 
pressure  distributions  exist  and  further  exper¬ 
imentation  and  exploration  are  needed  to  quantify 
the  nature  of  the  formation  of  film  rupture. 

6.  Careful  pressure  measurements  (from  moving 
and  stationary  surfaces)  of  both  the  statically 
and  the  dynamically  loaded  journal  bearing 
(although  it  is  a  difficult  task)  should  be 
complemented  by  film  temperature  measurements  in 
order  to  elucidate  some  of  the  unresolved  aspects 
of  cavitation. 
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CAVITATION  IN  JOURNAL  BEARINGS  AND  SQUEEZE  FILM  DAMPERS:  EXPERIMENTAL 
B.  JACOBSON2 


EXPERIMENTAL  INVESTIGATIONS 

In  nonstationary  bearings,  where  the  minimum 
film  thickness  moves  Jurlng  vibration,  vapor  cavi¬ 
tation  can  appear.  This  type  of  cavitation  Is  a 
breakage  of  the  oil  or  oil -bearing  interface 
caused  by  tensile  stresses  In  the  lubricant.  It 
can  appear  anywhere  in  the  oil  film  —  not  just  at 
the  outlet  where  the  film  thickness  is  Increasing 
In  the  direction  of  the  sliding  speed. 

Vapor  cavitation  has  recently  bet.i  studied 
by  Jacobson  and  Hamrock  (27,28)  and  Brewe  et  al . 
(29-31)  experimentally,  theoretically,  and  numeri¬ 
cally.  The  experiments  show  '•hat  the  cavltat’on 
Is  influenced  by  tne  vibration  amplitude,  the 
vibration  frequency,  and  the  sliding  velocity  In 
the  bearing.  The  clear  influence  of  the  surface 
energy  of  the  bearing  surfaces  Is  a’so  apparent; 
cavitation  seems  to  start  at  the  bearing  surface, 
If  It  Is  made  of  PTFE  (polytetrafluoroethylene) , 
or  in  the  middle  of  the  film,  If  the  surface  Is 
made  of  steel  or  PMMA  (polymethylmethacrylate) . 

In  experiments  vapor  cavitation  occurs  when 
the  tensile  stress  applied  to  the  oil  exceeds  the 
tensile  strength  of  the  oil  or  the  binding  of  the 
oil  to  the  surface.  The  physical  situation  neces¬ 
sary  for  vapor  cavitation  to  occur  Is  a  squeezing 
motion  within  a  bearing.  This  Implies  that  the 
minimum  film  thickness  Is  constantly  changing  In 
size  and  location  as  a  result  of  the  motion 
tangential  and  normal  to  the  surfaces. 

Experiments  showing  vapor  cavitation  were 
first  made  at  Lulea  Technical  Un.vers’-ty  In  the 
Machine  Elements  Laboratory  and  later  st  NASA 
Lewis  Research  Center. 

The  geometry  of  the  two  test  setups  was  the 
same,  but  the  NASA  Lewis  test  apparatus  was 
stlffer  and  more  precise,  both  In  the  measures  of 
the  bearings  and  also  In  the  motion  of  the  PMMA 
tube  and  the  shafts.  Figure  11  shows  the  NASA 
test  apparatus  with  the  high-speed  camera  and  the 
driving  motors  for  the  rotation  and  vibration 
motions.  The  high-speed  motion  picture  camera 
could  take  up  to  11  000  frames  per  second  with 
ordinary  16-mm  motion  picture  film.  In  the 
experiments  here,  the  framing  rates  were  1000  and 
2000  frames  per  second. 

The  test  bearing  was  lubricated  by  gravity 
feed  from  a  can  placed  approximately  600  mm  above 
the  bearing.  The  test  lubricant  was  a  nondegassed 
Dexron  II  automatic  transmission  fluid  with  a  vis¬ 
cosity  of  0.066  N-s/m2  under  atmospheric  condi¬ 
tions.  The  lubricant  passed  from  the  can,  through 


2Chalmers  University  of  Technology,  Goteborg, 
Sweden,  and  SKF  Engineering  and  Research  Centre, 
Nieuwegein,  The  Netherlands. 


Fig.  11 .  -  Test  apparatus. 


a  line  and  entered  the  bearing  at  the  top  left  of 
the  PMMA  bearing  (Fig.  11).  On  the  right  of  the 
test  bearing  is  the  motor  that  gives  the  sliding 
velocity  to  the  roller,  and  on  the  left  is  the 
motor  that  gives  the  squeeze  velocity  to  the 
bearing.  The  speeds  of  these  identical  0.5  kW 
motors  can  vary  continuously  from  0  to  250  rad/s. 

Figure  12  shows  the  PMMA  bearing,  which  can 
vibrate  In  a  circular  motion  keeping  the  axis  of 
vibration  parallel  to  the  axis  of  rotation.  A 
special  mechanism  was  manufactured  to  keep  the 
centerlines  of  the  bearing  and  shaft  parar->l  and 
Is  also  shown  in  this  figure. 

The  dynamic  eccentricity  mechanism  can  be 
adjusted  to  give  dynamic  eccentricity  ratios 
between  0  and  1  or  in  actual  measures  from  0  to 
0.5  mm.  The  static  eccentricity  ratio  is  set  by 
moving  the  motor  with  the  shaft  horizontally  so 
that  the  centerlines  of  the  two  motors  are  offset. 


Fig.  12.  -  PMMA-tube  and  eccentric. 
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Fig.  13.  -  Bearing  geometries  at  four  different 
times;  es  «  0.4;  ej  a  0.6 


(a)  (b)  (c)  (d)  (e)  (f)  (g) 

4  ms  12  ms  18  ms  20  ms  23  ms  24  ms  25  ms 


Fig.  14.  -  Dynamic  cavitation  bubble  developed 
(e$  *  0.34,  cd  =  0.66,  «s  »  19.5  r/s,  and  oitf  « 
92.7  r/s). 


vapor  cavitation  in  a  journal  bearing.  The  expert' 
mental  conditions  In  Fig.  14  are  length  to  diame¬ 
ter  ratio,  1/4;  static  eccentricity,  0.34;  dynamic 
eccentricity,  C.66;  rotatlonai  speed,  os,  19.5 
rad/s;  and  vibration  speed,  ud,  92.7  rad/s.  The 
journal  surface,  was  made  of  PTFE,  and  the  camera 
framing  rate  was  2000  frames  par  second.  The  time 
from  the  first  frame  to  the  last  frame  in  Fig.  14 
is  less  than  25  ms,  indicting  how  rapidly  the 
vapor  cavitation  forms  and  collapses. 


The  static  eccentricity  can  be  between  0  and  0.5 
mm.  The  offset  of  the  center  of  vibration  from 
the  center  of  the  journal  gives  a  static  load, 
and  the  amplitude  of  vibration  can  be  varied. 

Figure  13  shows  four  positions  of  the  journal 
relative  to  the  bearing  for  a  dynamically  loaded 
journal  bearing.  The  static  eccentricity  ratio 
es  is  set  at  0.4,  and  the  dynamic  eccentricity 
ratio  ed  Is  set  at  0.6. 

EXPERIMENTAL  RESULTS 

Figure  14  presents  seven  frames  from  the  high¬ 
speed  film  showing  the  formation  and  collapse  of 


Experiments  with  steel  journals  and  wl  h  PTFE 
journals  show  a  difference  in  the  cavitation 
pattern.  For  steel  journals  cavitation  bubbles 
start  In  different  parts  of  the  low  pressure 
region,  depending  on  where  the  largest  bubble  or 
dirt  particle  is  located.  For  PTFE  journals,  the 
cavitation  seems  to  start  from  the  same  position 
of  the  journal  for  each  vibration,  which  Indi¬ 
cates  a  rupture  of  the  oi 1 -PTFE  Interface. 

The  eccentricity  ratios  have  a  big  Influence 
on  the  size  of  the  cavitation  region  and  the  dura- 
tlo  of  this  region  from  start  to  collapse.  For 
small  dynamic  eccentricities  (compared  with  the 
minimum  oil  film  thickness),  no  vapor  cavitation 
Is  found. 
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CATEGORIZATION  OF  CAVITATION  TYPES  AND  ONSET  CONDITIONS 
D.H.  PARKINS3 


Oil  films  in  dynamically  loaded  journal 
bearings  are  subject  to  both  translatory  and 
oscillatory  squeeze  motion  of  their  containing 
surfaces.  This  discussion  deals  with  oscillatory 
squeeze  motion  only.  The  principal  aim  of  all 
work  described  here  is  to  understand  squeeze-film 
phenomena  -  simplified  by  removing  the  transla¬ 
tional  component.  It  is  important  to  distinguish 
sustained  oscillatory  squeeze  motion  from  "single¬ 
shot  movements"  in  which  the  oil  film  containing 
surfaces  approach  or  separate.  This  paper  deals 
with  the  former  case.  Some  published  work  <32,33) 
deals  with  oscillatory  motion  in  the  frequency 
range  of  1  to  2  Hz  and  only  the  transient  first 
few  cycles.  Other  works  consider  the  radial 
squeeze  motion  of  two  circular  surfaces,  as,  for 
example,  in  the  application  of  squeeze-film  dam¬ 
pers  to  the  outer  race  of  rolling-element  bear¬ 
ings.  This  discussion  deals  with  oscillatory 
normal  motion  of  parallel  surfaces  at  frequencies 
in  the  range  likely  to  be  met  in  machinery,  say  5 
to  50  Hz,  and  for  fully  established  behavior. 

Two  of  the  important  aspects  of  the  more 
recently  Identified  behavior  are 

1.  A  steady  force  is  generated  by  the  oil 
film  subject  to  normal  oscillatory  motion  of  its 
containing  surfaces.  The  possibility  was  first 
suggested  by  Kauzlarich  <34)  who  gave  a  theoret¬ 
ical  treatment  for  a  circular  pad  using  a  liquid 
lubricant  and  a  harmonic  displacement  input.  He 
concluded  that  inertial  forces  in  the  fluid  cre¬ 
ated  the  steady  force.  Kauzlarich  conducted  a 
qualitative  experiment  to  verify  this  finding. 
Subsequent  experimental  work  by  others  (35-37) 
using  nonharmonic  oscillatory  motion  of  a  square 
plate  has  confirmed  that  such  a  steady  force  is 
generated,  and  -  over  a  limited  range  of  param¬ 
eters  -  determined  the  Influence  thereon  of  the 
initial  oil  film  thickness,  the  oscillatory 
amplitude,  and  frequency  of  oscillation.  Most 
recently,  Kang  (38)  used  a  numerical  solution  ic 
the  Navier-Stokes  equation  to  compute  values  of 
the  steady  force. 

2.  Various  regimes  of  cavitation  appearing 
in  the  oil  film  due  to  the  oscillatory  squeeze 
motion  have  been  Identified  and  characterized, 
together  with  the  location  of  the  conditions  that 
determine  the  onset  of  cavitation.  Reference  39 
gives  an  initial  classification  and  Ref.  36  fuller 
descriptions,  all  for  behavior  of  an  oil  film  con¬ 
tained  in  a  square  plate.  Rodriques  (40)  predic¬ 
ted  the  growth  and  collapse  of  central  cavitation 
patterns  in  a  circular  plate.  His  supporting 
experiments  were  conducted  at  only  one  frequency, 
approximately  1.0  Hz  -  but  appear  to  show  cavita¬ 
tion  patterns  which  conform  exactly  to  the  charac¬ 
teristics  of  that  which  was  later  (39)  defined  as 
regime  2. 


3Cranfield  Institute  of  Technology,  College 
of  Manufacturing,  8edford,  MK43  Oal ,  England. 


STEADY  FORCE  GENERATED  BY  THE  OIL  FILM 

Figure  15  contains  data  on  the  magnitude  and 
direction  of  the  steady  force  <Fst)  for  eight 
cases  and  from  three  sources  (35,36,41).  Each 
source  used  the  same  basic  equipment,  but  other¬ 
wise  worked  completely  independently.  Some  curves 
show  data  obtained  with  invariant  vibratory  ampli¬ 
tude  (ha);  one  curve  shows  data  obtained  with 
invariant  frequency;  and  one  shows  amplitude  and 
frequency  varying  together.  The  data  of  Fig.  15 
are  for  initial  film  thicknesses  (y0)  from  0.11 
to  0.73  mm,  vibratory  amplitudes  from  0.0055  to 
0  5  mm,  frequencies  from  5  to  50  Hz,  three  values 
of  viscosity.  Two  values  of  oscillating  surface 
area  (A)  are  included.  The  use  of  the  steady 


Curve 

Ref. 

Amplitude 
of  oscillation 
ba. 
mm 

Initial 

film 

thickness, 

y0. 

mm 

Frequency, 

Pod 

area, 

mm2 

Lubricant 

kinematic 

viscosity, 

CS 

(a) 

35 

0.0055 

0.11 

10  to  50 

625 

313 

(b) 

35 

.0275 

.11 

10  to  50 

625 

313 

<c) 

41 

.044 

.22 

20  to  50 

2500 

750 

(d) 

41 

.044  to  .66 

.22 

10  to  50 

2500 

313 

(8) 

36 

.4 

.73 

5  to  45 

2500 

370 

(0 

36 

.4 

.5 

5  to  45 

2500 

370 

(g) 

36 

.4 

.25 

5  to  45 

2500 

370 

<h) 

36 

.07  to  .47 

.25 

25 

2500 

370 

Fig.  15.  -  Measured  steady  oil  film  pressure. 
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pressure  Pst  =  FS£/A  allows  for  its  effect. 
Reference  42  gives  additional  values  covering 
wider  ranges.  These  Include  data  for  a  third 
value  of  oscillatory  surface  area  A  and  from  a 
fourth  source  (43,44).  It  is  reassuring  to 
observe  that  all  the  curves  in  Fig.  15  form  a 
coherent  pattern,  even  though  they  originate  from 
separate  sources.  It  Is  also  reassuring  to  note 
that  the  additional  data  of  Ref.  42  would  also 
fit  the  coherent  pattern  if  Fig.  15  were  vastly 
extended.  However,  the  data  sparsely  cover  the 
parameter  range.  It  is  analogous  to  attempting 
to  determine  the  picture  of  a  jigsaw  puzzle  when 
only  a  few  scattered  parts  are  completed. 

Clearly,  there  remains  the  vital  task  of  con¬ 
ducting  a  complete  parametric  study  over  a  full 
range  of  variables. 

Data  from  Fig.  15  and  Ref.  42  show  that, 
depending  on  conditions,  the  steady  force  may 
either  increase  or  decrease  the  surface  separa¬ 
tion.  These  data  also  show  that  the  principal 
influences  on  the  steady  force  are  initial  film 
thickness  <y0),  oscillatory  amplitude  <ha),  and 
frequency.  The  steady  force  varies  directly 
with  ha  and  inversely  with  y0.  Cavitation 
appeared  in  many  of  the  cases  reported. 

CAVITATION  REGIMES 

It  is  widely  established  that  cavitation  may 
arise  in  such  oscillating  squeeze  films. 

Rodriques  (40)  measured  and  computed  cavitation 
patterns  with  a  circular  plate  oscillating  at 
1  Hz.  May-Miller  (39,41)  identified  three 
cavitatlng  behavioral  regimes,  which  he  named 
regimes  1,  2,  and  3.  Each  regime  is 
characterized  by  a  distinct  method  of  feeding  the 
bubbles  and  by  the  location  of  the  bubbles  in  the 
displacement  and  pressure  time  histories.  Exten¬ 
sion  of  the  concept  (36,45)  shows  that  regime  3 
could  take  several  forms.  Figure  16  shows  a 


FILM  THICKNESS 


Fig.  16.  -  Typical  oil  film  thickness  and 
pressure-time  history  (with  cavitation). 


typical  displacement  and  pressure  time  history 
with  cavitation  of  an  oscillatory  cycle  for  a 
harmonic  force  input.  The  pressure-time  history 
is  not  harmonic.  It  shows  a  rapid  increase  to 
maximum  pressure  followed  by  a  fall  to  the  mini¬ 
mum  value.  Regime  2  appears  in  this  part  of  the 
cycle.  A  very  rapid  rise  from  the  minimum  pres¬ 
sure  to  near  atmospheric  is  followed  by  a  rela¬ 
tively  long  period  near  atmospheric  pressure 
before  the  sharp  rise  again  to  maximum  pressure. 
Regime  3  appears  in  this  part  of  the  cycle.  The 
near  atmospheric  shoulder  in  the  pressure  curve 
of  Fig.  16  shows  a  small  positive  gradient  in  the 
versions  arising  when  regime  2  alone  appears. 

The  regimes  and  their  subdivisions  are 

Regime  1.  No  cavitation  at  any  time  in  the 
oscillatory  cycle. 

Regime  2.  Cavitation  bubbles  fed  from  within  the 
film.  Figure  17  is  a  photograph  of  a 
typical  fully  developed  regime  2 
cavitation  bubble  pattern.  The  bubbles 
may  contain  dissolved  gas  extracted 
from  the  solution  or  vapor  of  the 
liquid,  or  a  mixture  of  both. 

Regime  3.  Cavitation  bubbles  fed  from  outside  the 
film.  Cavitation  patterns  in  this 
regime  may  take  several  forms. 

(a)  Edge  fed:  This  form  is  characterized 
by  smoothly  curved  bubbles  entering 
simultaneously  from  all  edges  of  the 
square  plate.  A  typical  example  is 
shown  in  Fig.  18.  Some  bubbles  therein 
show  the  smooth  periphery  of  the  bubble 
breaking  into  strlations  which  lead  to 
the  form  described  in  under  3(b). 

(b>  Leaf  pattern:  This  form  is  character¬ 
ized  by  a  leaf  pattern  supplied  through 
a  "trunk"  that  grows  Inward  from  the 
plate  edge  (see  Fig.  19).  Usually  one, 


Fig.  17.  -  Regime  2  cavitation  pattern. 
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but  sometimes  more  such  trees  appear. 
Moreover,  the  tree  may  grow  from  a 
different  side  on  different  cycles  of 
the  same  condition.  References  36  and 
45  Illustrate  both  of  these  points. 

(c)  Leaf  pattern  with  square  parameter:  The 
leaf  pattern  extends  to  cover  all  the 
oscillating  plate  surface  except  for  a 
uniform  border  around  three  sides 

(see  Fig.  20). 

(d)  Stream  fed:  This  form  is  characterized 
by  a  leaf  pattern  whose  perimeter 
follows  that  of  the  oscillating  plate, 
but  with  a  relatively  wide  (b/4)  clear 
perimeter.  This  central  leaf  pattern  is 
fed  by  a  number  of  diagonal  streams. 
Figure  21  is  a  photograph  of  the  fully 
developed  pattern. 


Fig.  21.  -  Regime  3(d)  cavitation  bubble  pattern. 


Fig.  19.  -  Regime  3(b)  cavitation  bubble  pattern. 


Reference  36  gives  a  sequence  of  photographs 
of  cavitation  patterns  at  a  series  of  Identified 
points  in  the  oscillating  cycle  for  regimes  3(a), 
3(b),  and  3(d).  In  each  case  the  designated 
cavitation  pattern  exists  for  approximately  0.4  of 
the  cycle  time.  The  film  is  free  of  bubbles  for 
a  similar  period.  Bubbles  advance  and  withdraw 
in  the  remaining  time.  With  regime  3  (whatever 
subregime)  each  rapid  Increase  in  film  pressure 
arises  when  the  bubble  front  crosses  the  pressure 
transducer.  The  increase  from  minimum  to 
atmospheric  pressure  occurs  as  the  bubble  front 
advances,  and  that  from  atmospheric  to  maximum 
pressure  as  the  bubbles  withdraw.  An  important 
finding  in  this  work  (39,41)  was  the  extremely 
rapid  movement  of  the  bubble  front  in  both 
advance  and  withdrawal. 

Oil  depth  has  a  marked  effect  upon  the  onset 
cavitation  patterns  (36,45)  .  Figure  22,  which 
Illustrates  this  effect,  shows  the  value  of 
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Fig.  22.  -  Effect  of  oil  depth  on  onset  of  re¬ 
gime  3  cavitation  (Y0  a  0.7  mm). 


vibration  amplitude  (ha)  for  a  given  frequency 
and  oil  depth  (d)  at  which  regime  3  occurs.  Zones 
enclosed  by  dashed  lines  indicate  tvhlch  subregime 
will  appear.  For  example,  If  a  frequency,  say, 
of  20  Hz  Is  set  and  If  amplitude  Is  Increased 
from  zero,  then,  with  d  =  3  mm,  regime  3(a) 
would  appear  at  ha  =  0.57  mm,  regime  3(b)  at 
ha  =  0.75  mm,  and  regime  3(c)  at  ha  =  1.3  mm. 

It  is  Interesting  to  note  that  one  bubble  in 
Fig.  18  (regime  3(a))  shows  its  curved  edge 
changing  into  embryo  leaf  patterns  characteristic 
of  regime  3(b).  A  similar  exercise  at  a  frequency 
of  10  Hz,  d  =  3  mm  would  show  regime  3(a)  cavita¬ 
tion  at  ha  =  0.7  mm  with  no  change  therein 
occurring  as  amplitude  ha  is  increased 
further.  Alternatively,  for  a  frequency  of  20  Hz 
and  d  =  5  mm,  increasing  amplitude  ha  from 
zero  would  show  regime  3(b)  cavitation  at  onset 
at  ha  =  1 .2  mm. 

Figure  16  shows  how  regimes  2  and  3  occupy 
different  parts  of  the  cycle.  Consequently,  it 


is  possible  for  both  regimes  2  and  3  to  appear 
together  in  the  same  cycle.  Indeed,  Ref.  36 
gives  an  example  of  regimes  2  and  3(a)  appearing 
together.  Cavitation  onset,  of  whichever  regime, 
is  influenced  by  the  geometry  of  the  containing 
surfaces,  A,  y0,  and  vibration  frequency  and 
amplitude  ha.  However,  evidence  suggests  that 
the  onset  of  regime  3  is  markedly  affected  by  oil 
depth  (d).  The  onset  of  regime  2  however,  is 
influenced  by  the  dissolved  gas  content  in  the 
liquid.  These  are  independent  effects;  therefore, 
it  is  possible  to  have  many  different  combinations 
and  sequences.  It  is  possible  to  transform  from 
regime  1  (no  cavitation)  to  any  of  the  variants 
of  regime  3  with,  or  without,  regime  2  appearing. 
It  is  also  possible  to  transform  between  the 
subdivisions  of  regime  3.  Alternatively,  it  is 
possible  to  transform  from  regime  1  to  regime  2 
without  regime  3 ’appearing. 

This  large  number  of  possible  forms  of 
behavior  make  a  confused  picture  when  observed 
over  a  small  range  of  variables,  which  renders  it 
difficult  to  formulate  classifications.  There  is 
a  need  for  observation  to  be  taken  under  a  full 
range  of  all  variables. 


CONCLUDING  REMARKS 

It  is  now  well  established  that  an  oil  film 
subject  to  an  oscillatory  squeeze  motion  may  gen¬ 
erate  a  force  with  steady  and  vibrating  com¬ 
ponents.  The  steady  force  component  has  been 
measured  separately  and  shown  thereby  to  act  - 
according  to  circumstances  -  either  to  increase 
or  decrease  the  oil  film  thickness.  Theoretical 
models  that  include  fluid  inertia  have  been  shown 
to  predict  this  steady  force.  Measurements  have 
shown  that  the  initial  film  thickness,  vibratory 
amplitude,  and  frequency  influence  the  magnitude 
of  the  steady  force. 

Three  cavitation  regimes  -  including  sub¬ 
regimes  -  have  been  identified  and  characterized 
according  to  the  method  of  feeding  the  bubble. 

Each  regime  has  a  characteristic  form  of  displace¬ 
ment  and  pressure-time  history.  Regimes  2  and  3 
may  coexist  in  the  same  oscillatory  cycle.  A 
limited  knowledge  is  available  about  the  condi¬ 
tions  which  determine  the  onset  of  a  particular 
regime  or  subregime  and  its  subsequent  trans¬ 
formation  into  another  subregime. 

On  both  aspects,  steady  force  and  cavitation 
regimes,  current  knowledge  is  scanty,  but  suf¬ 
ficient  to  establish  the  existence  of  the  listed 
features.  More  work  needs  to  be  done  on  both 
fronts  before  a  full  picture  can  be  established. 
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CAVITATION  IN  MICROASPERITY  LUBRICATION 


A.O.  LEBECK* 


It  is  well  known  that,  even  when  sliders  are 
made  parallel  and  slide  parallel  in  a  lubricating 
film,  fluid  pressure  load  support  develops  so  as 
to  support  some  or  all  of  the  load  and  reduce 
friction.  Lebeck  (46,47)  reviewed  this  behavior, 
the  experimental  data,  and  the  mechanisms  In 
detai 1 . 

A  typical  example  of  such  behavior  is  shown 
in  Fig.  23.  The  friction  of  carbon  graphite  pins 
sliding  on  tungsten  carbide  In  water  was  measured 
using  the  apparatus  shown  In  Fig.  24.  The  fric¬ 
tion  coefficient  at  first  decreases  with  increas¬ 
ing  duty  parameter  and  then  later  Increases.  The 
decrease  Is  caused  by  what  are  apparently  hydro- 
dynamic  pressures  building  up  and  reducing  con¬ 
tact  friction.  The  later  increase  is,  of  course, 
caused  by  viscous  friction. 


Fig.  24.  -  Hydrostatic  bearing  assembly  diagram. 


For  many  years  it  has  been  thought  that 
microasperities  serve  as  a  source  of  fluid  pres¬ 
sure  load  support  in  mechanical  seals  and  parallel 
sliding  lubricated  contacts  as  reported  in  the 
experiments  by  Kojabashlan,  Richardson  (48),  and 
Hamilton  (49).  Figure  25  shows  how  mlcroasperl- 
tles  produce  load  support.  Cavitation  at  the 
downstream  side  is  essential  to  get  a  net  load 
support.  As  Kojabashlan  and  Richardson  (48)  show, 
one  can  predict  that  a  theoretically  significant 
load  support  will  develop  on  a  seal-like  surface 
which  has  many  microasperities. 

Experimental  tests  of  this  idea  show  that, 
if  the  ambient  pressure  is  raised  well  above  the 
cavitation  pressure,  no  net  load  support  results 
(Fig.  26).  It  follows  that  friction  should 


'////////////////////A 


'////////////, 


Fig.  23.  -  Parallel  sliding  in  water. 


••Mechanical  Seal  Technology  Inc.,  Albuquer¬ 
que,  NM. 


Fig.  25.  -  Truncated  asperity  pressure  distribu¬ 
tion. 


Fig.  26.  -  Zero  load  support. 
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Fig.  27.  -  Friction  torque  versus  ambient  pres¬ 
sure.  Carbon  pressure,  300  psi.  Each  point 
represents  an  average  of  450  samples. 


increase  because  a  lower  fraction  of  the  load  is 
being  supported  by  microasperities.  This  can  be 
verified  using  the  theory  of  Hamilton  et  al.  (50). 
To  test  this  hypothesis,  Kanas  (51)  ran  a  series 
of  experiments,  using  the  apparatus  shown  in 
Fig.  24,  where  friction  was  carefully  measured  as 
a  function  of  ambient  pressure.  His  results  (see 
Fig.  27)  show  no  significant  increase  in  friction 
with  increasing  ambient  pressure.  Similar  results 
were  obtained  at  higher  speeds. 

Thus,  in  water,  under  the  sliding  conditions 
noted,  and  at  a  condition  of  (from  Fig.  23)  con¬ 
siderable  fluid  pressure  load  support,  micro¬ 
asperity  lubrication,  or  any  other  mechanism  that 
relies  on  cavitation,  does  not  appear  to  be  a 
primary  source  of  load  support. 
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CAVITATION  EROSION  DAMAGE  IN  ENGINE  8EARINGS 


F.A.  MARTIN5 

The  experiences  discussed  in  the  1979  paper, 
"Cavitation  Erosion  Damage  in  Engine  Bearings: 
Theory  and  Practice,"  by  Garner,  James,  and 
Narriner  (52)  still  reflect  much  of  the  current 
view  on  the  subject.  A  few  selected  extracts 
from  that  fvaper  are  presented  here.  These 
include  examples  of  typical  cavitation  erosion 
damage,  together  with  suggested  mechanism  of 
damage.  Commonly  applied  palliatives  and  their 
effectiveness,  based  on  engine  experience,  are 
also  discussed. 

"Cavitation  erosion  in  engine  bearings  is  a 
phenomenon  which  has  assumed  increasing  signi¬ 
ficance  during  the  past. . .years ,  probably  as 
result  of  the  design  trends  towards  higher 
rotational  speeds  and,  in  some  cases,  higher 
rates  of  change  of  cylinder  pressure... 

"...Damage  is  seen  predominantly  in  diesel 
engine  bearings,  and  on  the  rare  occasions  when 
damage  occurs  in  gasoline  engines  It  is  normally 
as  a  result  of  operating  under  [incorrect 
conditions]. 

“In  many  cases  the  damage  is  restricted  to 
local  attack  of  the  overlay  and  the  performance 
of  the  bearing  Is  virtually  unaffected... 

However,  in  more  extreme  cases  extensive  loss  of 
overlay  and  interlayer  material  occurs,  and  the 
consequential  adverse  effects  upon  oil  film  cond¬ 
itions  will  reduce  the  service  reliability  and 
life  of  the  bearings. . . 

"The  fact  that  diesel  engine  bearings  are 
more  prone  to  cavitation  erosion  damage  Is  attri¬ 
buted  to  the  more  severe  combustion  conditions, 
and,  additionally,  to  the  necessity  for  more 
complex  oil  feeding  arrangements  with  greater 
chance  of  flow  discontinuities." 

Garner  et  al .  state  that  cavitation  erosion 
damage  is  associated  with  the  collapse  of  vaporous 
cavities  near  to  the  bearing  surface  and  suggest 
that  the  mechanism  is  the  impingement  by  micro¬ 
jets  of  fluid.  They  describe  a  possible  sequence 
of  events  in  the  collapse  of  vaporous  cavities 
leading  to  the  formation  of  high  intensity  micro- 
jets  (see  Fig.  28).  Let  us  now  look  at  the  result 
of  this  in  terms  of  actual  damaged  bearings. 

LARGE-END  8EARING  DAMAGE  PREDOMINANTLY  ASSOCIATED 
WITH  OIL  IN  CONNECTING  ROD  DRILLING 

"The  commonest  form  of  damage  occurring  in 
large  end  bearings  appears  at  first  sight  to  be 
due  more  to  fiow  erosion  than  to  cavity  collapse, 
but  detailed  examination  of  the  damage  confirms  a 
cavitation  effect.  This  can  occur  in  either 
fully  grooved  or  partially  grooved  bearings. 

'Tribology  Consultant  to  The  Glacier  Metal 
Company,  Limited,  Middlesex,  England. 
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Fig.  28.  -  Mechanism  of  micro-jet  formation  (52). 


"In  fully  grooved  bearings,  the  damage  occurs 
predominantly  at  the  groove  edges,  but  can 
ultimately  spread  into  the  bearing  lands.  There 
is  strong  evidence  that  the  damage  is  caused  by 
Inertia  effects  on  the  column  of  oil  in  the 
connecting  rod  drilling  which  is  supplied  from 
the  large-end  bearing  groove  to  provide  top  end 
lubrication  and  piston  cooling..." 

The  damage  may  be  cured  by  restricting  the 
back  flow  from  the  rod  drilling  into  the  bearing 
(fitting  a  nonreturn  valve  has  been  one  solution). 
In  other  cases,  "the  damage  was  diminished  when 
the  groove  cross-sectional  areas  in  both  the  main 
and  large-end  bearings  were  increased,  thus 
allowing  the  oil  column  fluctuation  to  be  better 
absorbed  by  the  oil  in  the  grooves.... 

"On  partially  grooved  large  end  bearings 
damage  frequently  occurs  at  the  groove  run-out, 
either  at  the  end  of  the  partial  groove,  or  at 
the  side  of  the  run-out,  [Fig.  29].  The  problem 
occurs  at  the  end  of  the  partial  groove  in  the 
direction  of  shaft  rotation,  regardless  of  the 
position  of  the  oil  feed  connection  to  the  rod 
drilling.  In  the  example  shown  in  [Fig.  29],  in 
which  the  damage  had  occurred  within  400  hours, 


Fig.  29.  -  Cavitation  damage  associated  with 
groove  run-out  in  large-end  bearing,  together 
with  effective  solution  (52). 
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the  oil  feed  to  the  rod  was  adjacent  to  the 
damaged  groove  ending.  The  problem  was  overcome 
by  the  simple  expedient  of  drilling  a  radial  hole 
in  the  run-out  section,  as  shown  by  the  lower 
bearing  in  [Fig.  29]  which  had  run  for  1500  hours. 
Several  other  cases  of  similar  damage  have  been 
successfully  treated,  ...either  by  provision  of  a 
radial  hole  or  by  terminating  the  groove 
conveniently  at  the  end  of  a  bridge-piece  slot." 


MAIN  BEARING  DAMAGE  PREDOMINANTLY  ASSOCIATED  WITH 
OIL  IN  CRANKSHAFT  DRILLING 


"The  majority  of  modern  medicm  speed  diesel 
engines  use  partially  grooved  main  bearings  in 
which  the  ungrooved  arc  is  symmetrically  placed 
in  the  bottom  half,  and  is  approximately  120®  in 
extent.  Such  bearings  provide  an  acceptable  oil 
film  thickness  within  the  narrow  confines  dictated 
by  the  overall  design,  but,  at  today's  high  ra¬ 
tings,  are  those  most  likely  to  suffer  from  cavi¬ 
tation  erosion  damage.  This  occurs  downstream  of 
but  usually  close  to,  the  groove  runout  and 
typically  takes  the  form  shown  in  [Fig.  30]. 

"In  this  instance,  the  damage  was  found  after  only 
400  hours  of  operation,  and  was  restricted  to  the 
soft  overlay  plated  surface  of  the  bearing,  which 
was  of  the  order  of  0.03  mm  thick.  The  stronger 
bearing  lining  underneath  had  not  been  attacked. 
This  form  of  damage  is  characterized  by  a  crescent 
form...  In  the  particular  example  shown,  after 
the  initial  short  term  damage,  further  progression 
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Fig.  31.  -  Cavity  formation  associated  with  inter¬ 
action  between  journal  drilling  and  partial 
groove  (52). 


was  slow.  Bearings  in  this  condition  after 
several  thousand  hours  operation  may  quite  safely 
be  refitted.  However,  in  some  Instances  the 
damage  can  be  much  more  severe... 

"It  is  probable  that  as  the  oil  drilling  in 
the  shaft  traverses  the  groove  run-out,  cavities 
are  formed  at  the  Intersection  of  groove  run-out 
edge  and  the  oil  hole  blend  radius,  as  depicted 
diagrammatical ly  In  [Fig.  31]." 

"Various  palliatives  have  been  tried  on 
different  engine  types  to  alleviate  this  type  of 
damage.  These  have  usually  consisted  of  modifi¬ 
cations  to  the  end  of  the  partial  groove.  The 
most  effective  palliative  so  far  adopted  has  been 
a  tangential  blend  from  groove  to  bearing  surface 
with  a  square  profile  ([Fig.  32])." 


Fig.  30.  -  Mild  overlay  damage  in  partially 
grooved  main  bearing  (52). 


Fig.  32.  -  Design  detail  of  successfully  applied 
modification  to  groove  ending  (52). 
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Fig.  33.  -  Suction  damage  In  large-end  bearing  as¬ 
sociated  with  large  clearance  <52). 


MAIN  AND  LARGE  END  BEARING  DAMAGE  PREDOMINANTLY 
ASSOCIATED  WITH  RAPID  SHAFT  MOTION 

Capitation  erosion  damage  of  the  type  shown 
In  C FI f  333  (commonly  referred  to  as  suction 
cavitat'on  erosion)  occurs  in  the  center  of  the 
lands  In  some  fully  grooved  bearings.  It  Is  due 
to  rapid  shaft  movement  away  from  the  bearing 
surface  and  Is  worse  with  large  clearances.  It 
occurs  in  the  upper  half  of  high-speed  diesel 
main  bearings  (and  also  In  a  large  end  bearing 
during  a  series  of  development  tests  on  a  medium 
speed  two  stroke  engine). 

Prediction  of  the  journal  orbit  and  exam¬ 
ination  of  the  hydrodynamic  conditions  within  the 
bearing  helps  to  pinpoint  the  area  of  likely 
damage.  It  also  enables  the  designer  to  quickly 
examine  the  relative  merits  and  disadvantages  of 
various  designs.  "For  the  purposes  of  balancing 
loads  and  integrated  pressures,  i.e.,  in  pre¬ 
dicting  the  journal  center  orbit,  it  has  been 
assumed,  as  Is  common  practice,  that  negative 
pressures  cannot  be  sustained  by  the  lubricant. 
However,  It  is  considered  that  predicted  negative 
pressures  wiii  indicate  the  position  and  inten¬ 
sity  of  any  cavity  formation." 

This  method  was  used  to  examine  the  hydro- 
dynamic  conditions  of  a  center  main  bearing 
showing  the  same  type  of  damage  as  shown  In 
Fig.  33.  "The  load  diagram  and  resulting  orbit 
are  shown  in  [Fig.  34(a)],  and  it  can  be  seen 
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Bearing  angle 
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Fig.  34.  -  Predicted  hydrodynamic  conditions  asso¬ 
ciated  with  occurrence  of  cavity  damage  (52). 


that  a  rapid  journal  movement  takes  place  between 
points  B  and  C.  Whilst  experience  suggests  that 
rapid  fluctuations  In  pressure  are  associated 
with  such  movements,  the  Initial  plot  of  rupture 
zones  within  the  bearing  ([Fig.  34(b)])  failed  to 
show  anything  significant  at  the  top  of  the 
bearing  (0  or  360°)  where  damage  occurred.  It 
was  only  when  the  time  step  of  the  solution  was 
refined  (from  10°  to  1°  crank  angle)  that  Intense 
"negative"  pressure  zones  were  predicted 
( CFfg.  34(c)])."  The  position  of  the  maximum 
negative  pressure  corresponded  to  the  position  of 
cavitation  damage  on  the  bearing  surface,  thus 
indicating  that  the  cause  of  the  damage  was 
cavity  formation  (and  subsequent  collapse).  The 
palliative  here  was  to  reduce  the  bearing  clear¬ 
ance  which  alleviated  the  damage.  However,  it 
may  be  necessary  in  other  Instances  to  modify  the 
polar  load  diagram.  For  example,  the  counter¬ 
weights  could  be  changed. 
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CONCLUDING  REMARKS 

Cavitation  erosion  damage  in  diesel  engine 
bearings  appears  to  fall  broadly  into  three  cate¬ 
gories,  associated  with 

1.  Oil  in  crankshaft  drillings  and  bearing 
oil  feed  features 

2.  Effects  of  the  oil  within  the  connecting 
rod  drillings 

3.  Journal  movements  within  the  bearing 
clearance  space  leading  to  adverse  oil  film 
conditions. 


Relatively  simple  modifications  to  the 
detail  design  of  the  bearing  and  Its  associated 
components  can  usually  obviate  problems  from  the 
first  two  categories,  but  often  not  from  the 
third.  However,  these  more  intractable  problems 
do  seem  to  show  some  correlation  with  predicted 
oil  film  conditions.  A  more  detailed  theory  than 
that  considered  here  will  probably  be  needed 
before  damage  can  be  reliably  predicted  at  the 
design  stage. 
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INTRODUCTION 

In  part  I,  we  were  made  aware  of  the  numer¬ 
ous  ways  that  cavitation  could  be  manifested  and_ 
of  the  categorization  thereof.  In  particular,  cav¬ 
itation  can  either  be  the  result  of  (a)  dissolved 
gas  coming  out  of  solution  or  (b)  evaporation 
(flashing)  of  the  fluid.  Both  types  of  cavitation 
are  commonly  observed  in  journal  bearings, 
squeeze-film  dampers,  connecting  rod  beamings,  and 
pistons.  Both  hav  a  pronounced  effect  ..,1  their 
operation.  The  oc  jrrence  of  cavitation  in  jour¬ 
nal  bearings  <1)  results  in  reduced  power  loss, 
friction  coefficient,  bearing  torque,  and  load 
capacity.  Dowson  and  Taylor  (2),  in  an  excellent 
review  of  cavitation,  point  out  that  cavitation 
need  not  have  a  deleterious  effect  on  the  load¬ 
carrying  capacity  of  bearings.  Horsnell's  calcula¬ 
tions  (3)  on  load  capacity  for  steady-state  condi¬ 
tions  show  that,  if  cavitation  Is  not  considered, 
the  load  can  increase  threefold  over  the  predicted 
load. 


In  recent  years  Imposing  higher  loads  and 
speed  together  with  more  complicated  loading 
cycles  resulted  in  bearings  experiencing  cavita¬ 
tion  erosion  damage,  as  discussed  in  the  article 
by  F.A.  Martin  in  part  I.  In  addition  to  the  cavi¬ 
tation  damage,  there  is  the  damage  that  can  occur 
due  to  self-excited  instabilities  that  are  encoun¬ 
tered  under  dynamic  loading.  These  instabilities 
can  be  manifested  as  a  whirling  or  whipping  motion 
(4)  of  the  journal  center.  Large  vibrational 
ampl'  udes  can  result  in  large  forces  being  trans¬ 
mitted  to  the  system.  A  designer  can  circumvent 
that  problem,  or  at  least  the  severity  of  It,  by 
knowing  the  speed  threshold  at  which  it  occurs. 

This  information  is  often  obtained  by  referring  to 
"stability  maps"  (5,6).  Stability  maps  require 
the  determination  of  hydrodynamic  force  terms  cou¬ 
pled  to  the  dynamical  equations  of  motion. 

The  calculation  of  the  hydrodynamic  force  com¬ 
ponents  is  dependent  on  the  film  model  used,  espe¬ 
cially  at  high  eccentricities.  Many  of  the  film 
models  that  are  used  to  generate  stability  maps 
involve  theories  that  are  overly  restrictive 
(e.g.,  narrow  bearing  theory,  infinite  bearing 
theory)  and  treat  the  cavitation  in  a  very  superfi¬ 
cial  way.  Notable  of  these  is  the  ir-film  cavita¬ 
tion  theory  in  which  positive  pressures  extend 
through  half  the  circumference  of  the  bearing  and 
the  other  half  is  regarded  as  cavitated.  This 
boundary  condition  violates  the  conservation  of 
mass.  Specifying  zero-pressure  gradient  and  cavi¬ 
tation  pressure  at  the  boundary  (i.e.,  Reynolds  or 
Swlft-Stieber  boundary  conditions)  represents  the 
conditions  for  film  rupture  in  a  reasonable  way 
but  does  not  properly  represent  the  conditions 
when  the  film  is  reestablished.  While  this  .ound- 
ary  condition  may  yield  satisfactory  result,  for  a 
steadily  loaded  bearing,  it  leads  to  inaccuracies 
in  situations  where  the  beat  ing  is  dynamically 
loaded  (7).  In  actual  practice,  the  Swift-Stieber 
conditions  are  often  compromised  as  a  matter  of 
convenience  by  neglecting  the  implementation  of 
the  zero-pressure  gradient  condition.  This  is 


much  1.  the  Giimbel  conditions  in  which  solutions 
were  ob  led  to  the  Reynolds  equation,  and  the 
cavitat  was  determined  by  disallowing  the  exist¬ 
ence  of  •  uambient  pressures.  However,  dynamic 
loading  uauses  changes  in  the  local  film  thick¬ 
ness,  which  leads  to  nonstationary  cavitation; 
that  Is,  the  cavitation  boundary  is  in  motion 
This  motion  is  manifested  as  the  growth  and 
collapse  of  the  bubble  as  well  as  downstream  trans¬ 
port  from  the  minimum  film  position.  The  appropri¬ 
ate  boundary  condition  may  require  a  condition 
other  than  a  zero-pressure  gradient  at  the  bound¬ 
ary.  The  studies  of  Olsson  (8,9)  suggest  that  the 
usual  Swift-Stieber  film-rupture  condition  is  ade¬ 
quate  for  dynamic  situations  if  the  cavitation 
boundary  moves  at  a  speed  that  is  less  than  half 
the  journal  surface  speed.  It  seems  highly  un”ke- 
ly  that  this  condition  would  be  met  during  the  Ini¬ 
tial  stages  of  growth  and  the  last  stages  of  bubble 
collapse. 


Jakobsson-Floberg  (1)  and  Olsson  (8) 
formulated  boundary  conditions  for  a  moving  bound¬ 
ary  that  conserved  mass  within  the  cavitated 
region  as  well  as  at  the  boundary  (commonly 
referred  to  as  the  JFO  cavitation  theory).  It  was 
assumed  that  liquid  was  convected  through  the  cavi¬ 
tated  region  in  the  form  of  strlations  extended  to 
both  surfaces  in  the  film  gap.  There  is  evidence 
that  some  liquid  is  transported  in  the  form  of  an 
adhered  layer  to  the  faster  surface  (10,11).  Pan 
(12)  has  broadened  the  JFO  theory  to  accommodate 
those  situations  in  which  the  mass  transport 
through  the  cavitated  region  is  not  necessarily 
via  liquid  strlations.  For  moderately  to  heavily 
loaded  bearings  (I.e.,  for  a  load-carrying 
capacity  much  greater  than  the  surface  tension 
forces  of  the  lubricant),  the  adhered  layer  can  be 
neglected  (10).  The  striated  flow  in  the  JFO 
theory  is  necessarily  a  Couette  flow  because  of  a 
constant-pressure  assumption  within  the  cavitated 
region.  Jakobsson  and  Floberg's  experimental  find¬ 
ings  (1)  support  the  constant-pressure  assumption. 


However,  Etsion  and  Ludwig  (11)  have  measured 
pressure  variations  of  the  order  of  50  kPa  inside 
a  gaseous  cavitation  bubble.  To  resolve  this 
apparent  contradiction,  one  must  know  whether  the 
cavitation  bubble  is  gaseous  or  vaporous.  Etsion 
(11)  provides  a  plausible  mechanism  for  determin¬ 
ing  this,  based  on  liberation  and  reabsorption 
rates  between  a  gas  and  a  liquid.  The  same  argu¬ 
ments  for  pressure  variation  within  a  vapor  bubble 
would  not  apply  (13).  In  view  of  this  discussion, 
the  constant  pressure  assumption  for  vapor  cavita¬ 
tion  seems  reasonable.  Despite  some  of  the  contro¬ 
versy  and  certain  lack  of  understanding,  the  JFO 
theory  perhaps  represents  one  of  the  best  accounts 
of  a  dynamical  theory  to  date  of  moderately  to 
heavily  loaded  journal  bearings  and/or  dampers. 

It  is  an  Improvement  over  the  Swift-Stieber  condi¬ 
tions,  even  for  steady-state  solutions,  because  it 
provides  for  film  reformation  as  well  as  film 
rupture.  Both  rupture  and  reformation  reouire  a 
knowledge  of  the  pressure  gradient  and  fractional 
film  content  at  the  interface  to  determine  its 
location.  Unlike  a  rupture  boundary,  the  reform- 
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ation  boundary  is  subjected  to  a  pressure  flow 
< t . e . ,  nonzero  pressure  gradient). 

Furthermore,  the  fractional  film  content  at 
the  boundary  is  determined  by  a  residual  fluid 
within  the  cavitated  region  that  has  been  released 
at  a  rupture  boundary  (earlier  in  time)  and  gov¬ 
erned  by  the  fluid  transport  law.  Which  condition 
(reformation  or  rupture)  prevails  depends  on  the 
relative  motion  of  the  boundary  with  the  motion 
of  the  convected  fluid  normal  to  the  boundary. 
Typically,  under  steady-state  cavitating  con¬ 
ditions,  the  rupture  boundary  occurs  along  the 
upstream  boundary,  and  the  reformation  boundary 
occurs  along  the  downstream  end.  Under  dynamic 
loading  and  nonstationary  cavitation,  rupture 
or  reformation  conditions  may  be  required  at 
either  the  upstream  or  the  downstream  boundary. 

The  programming  task  to  effect  these  conditions 
in  a  computational  algorithm  is  exceedingly 
tedious.  In  current  practice  this  discourages 
their  implementation. 

The  Elrod  and  Adams  (14)  computation  scheme, 
which  mimics  the  JFO  theory,  avoids  the  complex 
programming  required  to  trace  the  moving  boundary 


between  grid  points  and  to  evaluate  the  required 
pressure  derivative.  Later,  Elrod  (15)  modified 
it  and  presented  it  in  more  detail.  Rowe  and 
Chong  (16)  improvised  the  Elrod  algorithm  so  that 
pressure  Is  the  dependent  variable  rather  than 
fractional  film  content.  Provision  was  made  for 
pressure  sources  as  well  as  automatic  switching 
for  full  and  cavitated  regions.  Miranda  (17)  was 
one  of  the  first  to  apply  this  method  for  steadily 
loaded  axial  groove  journal  bearings.  The  full 
potential  of  the  algorithm  did  not  materialize 
until  Brewe  (7)  applied  It  to  nonstationary  cavi¬ 
tation  in  a  journal  bearing  undergoing  noncentered 
circular  whirl.  For  other  bearing  and  seal  appli¬ 
cations  using  the  Elrod  algorithm,  the  reader  is 
referred  to  in  (18-26) . 

This  work  represents  a  summary  of  the  state- 
of-the-art  and  is  by  no  means  complete.  Recently 
newer  ideas  and  methodologies  have  emerged  that 
contribute  to  the  present  methods  and  techniques 
of  numerical  modeling.  We  are  sincerely  grateful 
to  the  contributors  of  this  composition  in  present¬ 
ing  many  of  their  innovations,  ideas,  and  Insights 
to  this  complicated  problem. 
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NOMENCLATURE 


X]  First  Rlvl in-Erlcksen  tensor,  s-1 

A2  Second  Rlvl in-Ericksen  tensor,  s“2 

C  Bearing  radial  clearance,  n 

Cy  Damping  matrix,  N*s/m 

D  Bearing  diameter,  m 

e  Eccentricity,  m 

F  Total  fluid  force,  N 

F(s)  Deformation  gradient  history 

f  Dynamic  fluid  force,  N 

f/\  Adhered-fllm  fill  factor 

g  Switch  function 

H  Nondlmens lonal  film  thickness,  (h/C) 

Ha  Nondlmenslonal  adhered  film  thickness 

h  Film  thickness,  m 

he  Film  coeflclent  of  heat  transfer  for  condens¬ 
ing  oil  vapor  at  low  flow  velocity,  W/m*  C2 

I  Last  grid  point 

1  Grid  point  Index 

Kij  Stiffness  matrix,  N/m 

L  Length  of  bearing  In  z-dlrectlon,  m 

L'  Axial  halflength  (L1  ■  1/2L),  m 

lv  Heat  of  vaporization,  kJ/kg 

N  Total  number  of  time  steps 

n  Iteration  level 

n  Unit  normal  vector  of  the  Internal  boundary 

Ob  Bearing  center 

Oj  Journal  center 

P  Nondimenslonal  pressure,  p/pu(C/R)2 

p  Fluid  pressure,  N/m2 

Pvap  Vapor  pressure,  N/m2 

pa  Atmospheric  pressure,  N/m2 

6  Sweeping  flux  vector,  m2/s 

QZL  Nondlmenslonal  axial  flow 


R  Bearing  or  journal  radius,  m 
r  Radius  of  flat  bubble,  m 

s  Marking  parameter  In  time,  s 

T  Temperature,  °C 

ATvj  Temperature  change  from  oil  vapor  to 
liquid,  °C 
T  Stress  tensor 

t  Time,  s 

U  Surface  velocity  In  x  direction,  and/or 
sliding  velocity,  m/s 
u  Velocity  of  the  fluid  In  the 

circumferential  direction,  m/s 

V  Collapse  velocity  of  the  bubble,  m/s 

v  Velocity  vector,  m/s 

W  Applied  static  load  on  bearing,  N;  width  of 

cavitation  bubble,  m 

X  Nondlmensional  x-coordlnate,  x/R 

x  Coordinate  In  circumferential  direction,  m; 

Instantaneous  position  coordinate  relative 
to  journal  center,  m 

xs  Vertical  displacement  of  journal  center  rel¬ 
ative  to  bearing  center,  m 

x  Instantaneous  position  coordinate  relative 

to  bearing  center,  m 

x  Instantaneous  x-component  of  journal  veloc¬ 

ity,  m/s 

Y  Nondlmenslonal  y-coordinate,  y/R 

y  Coordinate  In  radial  or  normal  direction,  m; 
Instantaneous  position  coordinate  relative 
to  journal  center,  m 

y  Instantaneous  y-component  of  journal  veloc¬ 

ity,  m/s 

y  Instantaneous  position  coordinate  relative 

to  bearing  center,  m 

ys  Horizontal  displacement  of  journal  center 
relative  to  bearing  center,  m 

Z  Nondlmenslonal  z-coordlnate,  z/L * 

z  Coordinate  In  axial  direction,  m 
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zc  Location  of  bubble,  m 

ze  Axial  location  of  the  cavitation  boundary,  m 

zm  1/2  Bearing  width,  m 

a  Inclination  of  unit  normal  (n)  to  direction 

of  sliding,  rad;  non-Newtonian  material 
parameter,  N-s/m2 

(J  L'quid  bulk  modulus,  N/m2;  non-Newtonian 

material  parameter,  N*s2/m2 

Y  Nondimensioned  counterpart  to  the  material 

parameter  a,  (aU/pR) 

S  Central  difference  operator;  nondimensioned 

counterpart  to  the  material  parameter  (J, 
(6|3U2/pC2) 

c  Eccc’ii’ricity  ratio  of  bearing  (e/C) 

5  Position  of  upstream  cavitation  boundary  at 

any  axial  position,  m 

n  Unwrapped  circumferential  coordinate  along 

the  cylindrical  bearing  surface,  m 

0  Nondimensional  density  in  full  film,  <p/pc); 
fractional  film  content  in  cavitation  region 

0n  Fill  factor  of  adhered  layer  on  sliding 
surface 

9(3  Bridged  fill  factor 

Qq  Fill  factor  of  layer  clung  to  stationary 
surface 

X  Ratio  of  bearing  half  length  to  bearing  ra¬ 
dius,  <L7R0> 

p  Viscosity,  N-s/m2 

pa  Artificial  viscosity,  N*s/m2 

5  Positional  variable,  m;  Axial  coordinate  of 

the  cylindrical  bearing  surface,  m 

P  Density  of  saturated  oil  vapor,  kg/m3 


t  Period  of  dynamic  load;  the  instant  the  ad¬ 
hered  film  was  formed  at  the  breakup  point; 
former  instant  in  time,  s 

xXy  Shear  stress,  N/m2 

Y  Attitude  angle,  rad. 

<i>  Film  flux,  m2/s 

X  Inverse  function  for  the  spatial  position  x, 
m 

Xt  Relative  deformation  function 

$  Liquid  phase  transport  flux,  m2/s 

Ye  Circumferential  transport  flow,  m3/s 

Yz  Axial  pressure  flow,  m3/s 

«  Angular  velocity  of  shaft,  rad/s 

Subscripts 
a  Surface  a 

b  Surface  h 

B  Bearing 

c  Cavitation 

d  Dynamic 

F  Full  fluid  film 

j  Journal 

max  Maximum 

min  Minimum 

S  Starvation 

s  Static,  supply 

w  Waist 
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DEVELOPMENT  AND  EVALUATION  OF  A  CAVITATION  ALGORITHM 


T.G.  KEITH'  AND  D.  VIJAYARAGHAVAN' 

The  effects  of  cavitation  on  performance  have 
largely  been  disregarded  In  bearing  performance 
calculations.  For  years,  common  practice  was  to 
modify  full-film  results  in  load  calculations  by 
setting  all  predicted  negative  pressures  to  the 
cavitation  pressure  (the  half-Sommerfeld  or  GUmbel 
boundary  conditions).  Albeit  this  produced  reason¬ 
able  load  values,  it  is  also  well  known  that  the 
results  violated  the  conservation  of  mass  princi¬ 
ple.  In  more  recent  procedures  the  film  pressure 
is  assumed  to  develop  from  the  point  of  maximum 
film  thickness  and  end  where  the  pressure  deriva¬ 
tive  vanisnes  (the  Reynolds  or  Swift-Stieber  bound¬ 
ary  condition).  Film  rupture  is  more  appropriately 
treated,  but  reformation  of  the  film  Is  essentially 
neglected. 

Jacobsson  and  Floberg  (1)  and  Olsson  (7),  pro¬ 
posed  a  set  of  self-consistent  boundary  conditions 
for  cavitation  to  be  applied  to  the  Reynolds  equa¬ 
tion.  The  boundary  conditions  properly  account 
for  mass  conservation  in  the  cavitation  region. 

The  collective  procedure  is  now  generally  termed 
JFO  theory.  This  theory,  when  applied  to  moder¬ 
ately  to  heavily  loaded  journal  bearings,  produced 
values  that  were  in  satisfactory  agreement  with 
experimental  values.  The  theory  markedly  improved 
the  prediction  of  cavitation  regions  in  bearing 
analyses  and  for  transient  problems. 

Elrod  and  Adams  (13)  introduced  a  computa¬ 
tional  scheme  that  Incorporates  JFO  theory  In  a 
very  simple  manner.  This  procedure  avoids  the  com¬ 
plexities  of  locating  film  rupture  and  reformation 
boundaries  as  It  automatically  predicts  cavitation 
regions.  Later  Elrod  (14)  modified  the  finite  dif¬ 
ference  portion  of  the  procedure  to  provide  better 
results.  The  procedure,  now  generally  termed  the 
Elrod  cavitation  algorithm,  has  the  following  five 
objectives: 

1.  To  remove  the  pressure  contribution  to  the 
mass  flow  In  the  cavitation  region 

2.  To  observe  the  conservation  of  mass 

3.  To  set  approximately  correct  boundary 
conditions  at  the  boundaries  of  the  cavitation 
region 

4.  To  produce  a  single  equation  for  both  the 
full  film  and  the  cavitated  regions 

5.  To  Introduce  the  proper  finite  differenc¬ 
ing  of  the  governing  equation. 

The  algorithm  is  based  on  a  control  volume 
formulation  and  uses  the  unit  step  (or  switch) 
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function  to  eliminate  the  pressure  term  in  the 
cavitation  region.  Elrod  pointed  out  that  the 
algorithm  was  developed  "after  considerable 
experimentation"  and  therefore  does  not  present 
development  details.  Accordingly,  it  does  not 
readily  lend  itself  to  improvement. 

An  attempt  is  made  here  to  develop  an  algo¬ 
rithm  that  Is  as  effective  as  the  Elrod  algorithm 
in  predicting  cavitation  In  bearings.  The  pro¬ 
posed  algorithm  has  the  same  objectives  as  the 
Elrod  algorithm  but  Is  developed  without  resorting 
to  trial  and  error.  Hopefully,  this  development 
will  serve  as  a  platform  upon  which  further  algo¬ 
rithm  improvement  can  occur. 

ANALYTICAL  FORMULATION 

The  two-dimensional,  transient  form  of 
Reynolds  equation  for  a  Newtonian  lubricant  in  lam¬ 
inar  flow,  allowing  for  compressibility  effects, 
can  be  written  as 


This  is  an  expression  of  the  conservation  of 
mass:  the  first  term  on  the  left  is  the  squeeze 
term;  the  second  is  the  net  mass  flow  rate  In  the 
x  (flow)  direction;  and  the  third  is  the  net  mass 
flow  rate  In  the  z  (axial)  direction.  Note  that 
the  flow  rate  in  the  x  direction  consists  of  a  con¬ 
tribution  due  to  shear  (Couette  flow)  and  a  contri¬ 
bution  due  to  pressure  gradient  (Polseullle  flow); 
whereas,  the  flow  rate  In  the  z-dlrectlon  is 
solely  due  to  pressure. 

The  density  of  the  lubricant  is  related  to 
the  film  pressure  through  the  definition  of  the 
bulk  modulus,  P,  i .e. , 

P  =  p(3p/3p)  C2] 

This  relation  suggests  that  It  Is  easier  to  work 
in  terms  of  density  in  Eq.  [1]  than  in  terms  of 
pressure.  Accordingly,  we  can  define  a  non- 
dimensional  density  variable 

0  =  p  =  pc  [3] 

where  pc  is  the  density  of  the  lubricant  at  the 
cavitation  pressure. 

For  moderate  to  heavily  loaded  bearings,  the 
pressure-induced  flow  within  cavitated  regions  is 
small.  Thus,  the  striated  flow  existing  in  these 
regions  is  primarily  driven  by  shear,  and  the  pres¬ 
sure  is  essentially  constant.  To  reflect  this  fea¬ 
ture,  a  so-called  switch  function  is  Introduced 
into  the  pressure-density  relation  of  Eq.  [2], 
l.e., 

gP  =  p(3p/3p)  =  9(3p/3p)  [4] 
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where 

II  in  the  ful  1  fi  lm  region 

0  in  the  cavitasea  region 

Equation  [4]  can  be  d,.ectly  integrated  to  yield 

p  =  pc  +  g(i<  1  n  6)  [5] 

Using  Eqs.  Cl],  C 3 3 ,  and  C4]  allows  use  to  write  a 
single  expression  for  the  conservation  of  mass  in 
the  bearing 


In  order  to  reflect  the  correct  physics  of 
the  flow,  we  desire  a  central  difference  of  this 
term  In  the  full-film  portion  of  the  bearing  and  a 
one-sided  backwards  difference  in  the  cavitated 
region.  This  can  be  accomplished  automatically  by 
the  insertion  of  an  artificial  viscosity  term  into 
Eq.  C8].  This  procedure  is  called  type  differenc¬ 
ing  (27)  when  used  for  similar  purposes  in  tran¬ 
sonic  flow  computations.  Thus,  Eq.  [8]  is  written 
as 


3pche  3  /pchu  pcph39  ae 
at  +  ax  [  2  0  "  i2p  ax 
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In  the  cavitated  region  Eqs.  <5)  and  (6)  are 

8.  (P  hue)  =  0  p  =  p 
at  ax  u  ;  c 


[63 
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The  variable  6  warrants  some  discussion, 
for,  although  it  is  defined  as  the  ratio  of  densi¬ 
ties,  it  carries  a  somewhat  different  interpreta¬ 
tion  In  the  full  film  than  in  the  cavitation 
region.  In  the  former  e  >  1  and  p  >  pc  because 
of  the  compression  of  the  lubricant.  In  cavitation 
regions  0  <  1  for  another  reason:  The  density  of 
the  lubricant  is  uniformly  pc,  but,  because  of 
the  vaporization,  the  film  Is  not  completely  occu¬ 
pied  with  lubricant.  Thus,  0  Is  termed  the  frac¬ 
tional  film  content  and  1-0  represents  the 
void  fraction. 


where  pa,  the  artificial  viscosity  switch  func¬ 
tion,  is  defined  as 


0  in  the  full  film  region 
1  in  the  cavitated  region 


The  use  of  the  switch  function  permits 
Eq.  C93  to  be  centrally  differenced  to  yield 
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Values  of  the  variables  on  the  half  step  may  be 
computed  from 


NUMERICAL  FORMULATION 


Examination  of  the  steady-state  version  of 
Eq.  [63  for  the  full  film  (g  =  1)  reveals  that  it 
is  an  elliptic  partial-differential  equation. 
Hence,  the  numerical  scheme  adopted  for  its  solu¬ 
tion  should  reflect  the  fact  that  the  dependent 
variable  at  any  point  in  the  full  film  depends  on 
all  neighboring  variables.  This  permits  the  use 
of  central  differences.  On  the  other  hand,  in  the 
cavitated  region  <g  =  0),  Eq.  [63  is  of  the  hyper¬ 
bolic  type.  Accordingly,  downstream  influences 
are  not  signalled  upstream,  and  thus  one-sided  or 
upwind  differencing  should  be  used. 

In  the  following  discussion,  the  shear  flow 
and  the  pressure-induced  flow  components  will  be 
treated  separately.  A  one-dimensional  flow  will 
be  used  in  order  not  to  obscure  the  details  of  the 
development.  Extension  to  multidimensional  flow 
is  not  difficult. 

Shear  Flow 


For  simplicity  the  shear  flow  term  will  be 
written  as 


(  ) 
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i+1 72 


Using  Eq.  [11],  performing  central  differencing  on 
the  remaining  derivatives  in  Eq.  [lu],  and 
rearranging  produce 


6_  ,U  aEi  +  l  *  bEl  *  cE1-1  [12] 

Sx  (2  pch  l  =  2  Ax 


where 


a  =  1  -  (pa>i+l/2 
b  =  Wl+1/2  +  Wi-1/2 
c  =  -[1  +  Wl-1/23 


[13] 


Comparison  of  the  definition  of  pa  and  the 
switch  function  introduced  to  cancel  the  pressure 
term  in  the  cavitated  region,  via  Eq.  [4],  reveals 
that 


g  =  1  -  Pa  [14] 

Thus,  the  coefficient  in  Eq.  [133  are 
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a  =  (gi+i  +  gj)/2 


b  =  2  -  (g)  +  ]  +2gi  +  gj_])/2  >•  C153 


c  =  -[2  -  (gj  +  gi_])/2] J 

Now,  If  all  grid  points  are  within  the  full  film, 
all  g's  are  unity;  therefore,  a  =  1 ,  l  =0,  and 
c  =  -1 .  Hence, 


This  term  exists  in  the  full  film  region  (6  >  1), 
but  vanishes  In  cavltated  regions  <9  <  1).  The 
switch  function  g  adequately  removes  the  term 
and,  therefore,  no  artificial  viscosity  term  need 
be  added  for  this  purpose.  Also,  when  this  term 
exists,  the  governing  equation  Is  elliptic  and 
thus  can  be  centrally  differenced. 


pcu  f<h©)i+1  -  (h0),_, 
__ 
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which  Is  the  desired  central  difference  form.  On 
the  other  hand,  If  all  grid  points  are  within  the 
cavltated  region,  all  g's  vanish  so  that  a  =  0, 
b  =  2  and  c  =  -2.  Hence, 


Equation  [213  can  be  treated  In  a  variety  of 
ways  which  will  ensure  that  it  Is  central  differ¬ 
enced  In  the  full-film  and  made  to  vanish  In  the 
cavltated  regions.  However,  depending  on  the  man¬ 
ner  In  which  the  switch  function  Is  handled  and 
the  particular  differencing  scheme  adopted,  the 
results  vary  significantly  at  the  cavitation  bound¬ 
aries.  Several  alternative  versions  to  the  above 
were  examined.  The  version  that  provided  the  most 
adequate  and  consistent  results  Is  one  In  which 
the  derivative  of  g  Is  not  explicitly  differ¬ 
enced.  In  particular,  we  may  write  that 


E,  -  E, 
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and  the  desired  upwind  difference  version  Is 
obtained. 

The  foregoing  may  be  contrasted  to  Elrod's  (15) 
shear  flow  expression,  which  Is  written  as 

fx(l  pch6),  3  "  2Sx  hi+l(1  “  gi-l)01-l 
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If  all  grid  points  are  within  the  full  film,  this 
yields 


.30  .3(0-1)  3q(9  -  1 )  , ,3q 

g__=  g^_ -  =  - (0  "  ^ax  C22] 

The  last  term  on  the  right  hand  side  vanishes 
for  all  values  of  9  because  the  spatial  deriva¬ 
tive  of  g  vanishes  for  all  9  except  at  9=1; 
however,  at  this  point,  obviously,  9-1  Is 
zero.  Hence,  we  have 

39  3g  (9  -  1)  [233 


Now  Inserting  Eq.  [233  Into  Eq.  [21]  and  central 
differencing  yields 
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which,  curiously,  does  not  account  for  compressi¬ 
bility  effects.  In  the  cavitation  region, 

Eq.  [183  produces  r  , 
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which  agrees  with  Eq.  [173. 


Pressure  Induced  Flow 

The  pressure  Induced  flow,  written  in  terms 
of  9  Is 


x  ^*1-1/2  gi-l^0,-l  "  £24] 

which  Is  the  same  pressure  term  expression  used  by 
Elrod  (15). 


SOLUTION  PROCEDURE 

In  this  study  cavitation  In  both  slider  and 
journal  bearings  was  examined.  All  sliders  were 
assumed  to  be  one  dimensional,  whereas  both  finite 
and  infinite  length  journal  bearings  were  consid¬ 
ered.  Because  the  procedures  for  solving  one-  and 
two-  dimensional  problems  in  these  configurations 
differed,  details  of  each  will  be  briefly  described 
In  the  following  sections. 
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One-Dimensional  Problems 


A  purely  implicit  finite-difference  scheme 
was  used  to  numerically  determine  the  distribution 
of  9  for  all  one-dimensional  problems.  By  this 
method,  the  finite  difference  version  of  Eq.  [6], 
which  Incorporates  Eqs.  [12],  [15],  and  [24],  can 
be  written  at  each  node  in  a  general  form  as 


a  fln+' 

aiei-1 
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where  the  coefficients  a),  bj ,  cj ,  and  d^  are 
functions  of  Ax,  At,  |5,  p,  h,  g,  and  the  form  of 
the  boundary  conditions  applied.  Consequently, 
the  coefficients  vary  throughout  the  grid  and 
change  with  time.  Initially,  all  nodes  are  as¬ 
sumed  to  be  in  the  full  film,  and  at  all  nonboun¬ 
dary  nodes  6  =  1  and  g  =  1.  The  switch  function 
distribution  was  updated  after  every  time  step. 

Boundary  conditions  for  slider  bearings  were 
taken  to  be  of  the  Dirichlet  type  in  which  9  was 
set  either  greater  than  one  (flooded)  or  less  than 
one  (starved)  on  the  boundary  nodes.  Because  the 
inlet  and  exit  values  of  9  were  specified,  the 
entire  collection  of  nodal  equations  could  be 
arranged  into  the  following  tridiagonal  matrix 
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be  directly  used  in  the  solution.  Brewe  (7)  has 
used  the  Gauss-Jordan  procedure  for  this  purpose, 
but  found  it  to  be  extremely  slow.  To  circumvent 
this  problem,  the  two  off-diagonal  terms  were 
lagged  during  the  computations.  By  this  procedure 
the  following  replacements  were  made: 
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Thus,  the  terms  were  evaluated  one  time  step  be 
hind  the  actual  computations  and,  therefore,  are 
known  during  the  (n+l)st  time  step.  Accordingly, 
they  were  Incorporated  into  the  right  hand  side  of 
Eq.  [27],  and  the  Thomas  algorithm  applied. 


The  time  march  process  was  continued  until 
the  change  In  9]  between  two  time  intervals  was 
less  than  10-4.  At  this  point,  it  was  assumed 
that  a  converged  steady  state  had  been  reached. 


Two-Dimensional  Problems 


An  alternating  direction  implicit  (ADI) 
finite  difference  scheme  (7)  was  used  to  determine 
the  distribution  of  9  for  all  finite  journal 
bearing  problems.  By  this  method,  the  time  step 
is  split  into  two  parts:  in  the  first  half, 
changes  in  the  flow  or  x  direction  are  calculated 
while  old  9  values  are  used  to  determine  the  gra¬ 
dients  in  the  axial  or  z  direction;  in  the  second 
half,  the  recently  calculated  9  distribution  is 
used  to  evaluate  the  flow  direction  gradients 
while  the  axial  changes  are  updated. 


[20] 

In  this  form  the  collection  of  equations  may  be 
readily  solved  by  application  of  the  Thomas 
algorithm  (27). 


For  journal  bearings  the  boundary  conditions 
depend  on  the  means  of  lubricant  supply.  At  oil 
grooves,  values  of  9  are  assumed  which  corre¬ 
spond  to  the  given  supply  pressure  and  are  calcula¬ 
ted  from  Eq.  [5].  The  procedure  remains  the  same 
as  for  slider  bearings.  If,  however,  the  bearing 
is  submerged,  periodic  or  continuous  boundary  con¬ 
ditions  must  be  employed.  Application  of  these 
conditions  produce  the  following  matrix  form  of 
nodal  equations 
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The  coefficient  matrix  in  Eq.  [27]  is  not  tri¬ 
diagonal  and  therefore  the  Thomas  algorithm  cannot 


The  Thomas  algorithm  can  be  used  in  each  part 
of  the  time  step,  If  again,  the  appropriate  varia¬ 
bles  are  lagged  for  the  submerged  journal  bearings 
case.  In  accordance  with  Brewe1 s  (7)  suggestion, 
the  switch  function  was  updated  after  every  half 
time  step.  Calculations  were  performed  over  half 
the  length  of  bearing  due  to  symmetry  with  resepct 
to  the  axial  center. 

On  occasion  an  oscillation  of  the  9  distri¬ 
bution  was  encountered.  This  oscillation  was  trig¬ 
gered  at  the  cavitation  node  and  has  been  reported 
elsewhere  (e.g.,  (7)).  To  remedy  this  problem,  a 
variable  time  step  was  implemented  according  to 
the  following  scheme: 


■  atl,n/at»“Vk~'>,<IM> 


min, 


k  =  1,  2 .  N 

[29] 


where  Atmin  and  Atmax  are  grid  dependent  and 
N  is  normally  taken  to  be  8. 


Not  only  did  the  variable  time  step  help 
reduce  the  oscillation,  it  also  reduced  by  approxi¬ 
mately  15  percent  the  number  of  time  steps  re¬ 
quired  for  the  solution  to  reach  steady  state. 
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p  *  pc  +  (i(©  -  1) 


Reference  28  examines  the  details  of  the  cur¬ 
rent  approach  for  both  slider  and  journal  bearings 
in  one  and  two  dimensions.  The  principal  purpose 
of  that  effort  was  to  compare  the  computed  results 
with  the  published  values  of  Elrod  (14,15):  The 
two  algotithms  produced  comparable  results.  A  sec¬ 
ond  purpose  of  that  work  was  to  examine  a  few  modi¬ 
fications  to  the  basic  flow  terms  used  in  the  Elrod 
algorithm.  In  particular,  Elrod  neglects  in  the 
shear  flow  term  the  compressibility  effects  in  the 
full  film  region.  Also,  an  approximate  equation 
relating  the  pressure  to  the  fractional  film  con¬ 
tent  is  sometimes  used.  This  approximate  relation 
uses  the  first  two  terms  of  a  series  expansion  of 
the  logarithmic  term  In  Eq.  C5],  i . e . , 


[30] 

CONCLUSIONS 

A  cavitation  algorithm  that  automatically 
alters  the  form  of  finite  differencing  used  for 
the  shear  Induced  flow  has  been  developed.  The 
algorithm  is  an  outgrowth  of  the  Elrod  cavitation 
algorithm.  The  two  algorithms  produce  comparable 
results  (28).  However,  the  newer  version  may  lend 
Itself  more  readily  to  Improvement.  The  effects 
of  film  compressibility,  equation  of  state  of  the 
lubricant,  and  grid  spacing  near  cavitation  bounda¬ 
ries  on  the  pressure  distribution  were  investi¬ 
gated  (28). 


33 


FUNDAMENTAL  ISSUES  CONCERNING  THE  MOVING  BOUNDARY 


CODA  H.T.  PAN2  AND  JEAN  FRENE3 

This  discussion  deals  with  the  motion  of  an 
internal  cavitation  boundary  of  an  Incomplete 
film,  with  emphasis  on  applications  in  which  the 
dynamic  load  component  overshadows  the  magnitude 
of  the  static  load.  In  such  applications  one  is 
concerned  with  supplying  lubricant  to  the  region 
where  the  surfaces  are  approaching  each  other, 
while  voids  should  dominate  the  region  of  separat¬ 
ing  surfaces  such  that  viscous  friction  may  be 
kept  minimal.  The  commonly  accepted  model  of  the 
tt  film  is  an  attempt  to  approximate  the  extent  of 
filling  of  the  bearing  gap.  As  engineers  desire 
to  refine  the  design  art,  a  more  precise  descrip¬ 
tion  of  the  internal  boundary  of  an  Incomplete 
film  In  a  bearing  gap  Is  sought. 

THE  MODELING  VIEWPOINT 


A  predominant  concern  In  the  analysis  of  an 
Incomplete  film  is  related  to  how  flow  continuity 
is  maintained  across  the  entire  perimeter  of  the 
film,  which  comprises: 

1.  Inlet  film  edge  or  edges. 

2.  Boundaries  fixed  by  the  geometry  of  the 
bearing. 


3.  Internal  boundary  or  boundrles. 

Suppose,  that  at  any  Instant,  such  delinea¬ 
tion  of  the  film  perimeter  can  be  asserted.  Then 
one  can  assign  the  value  of  the  film  pressure  at 
the  film  Inlet  to  be  that  of  the  lubricant  supply 
pressure;  the  geometrical  boundaries  can  be  repre¬ 
sented  by  the  ambient  pressure;  and  the  "cavity" 
pressure  can  be  assumed  to  prevail  near  the  Inter¬ 
nal  boundary  or  boundrles.  Consistent  with  this 
set  of  pressure  boundary  conditions,  a  Dlrichlet 
type  problem  Is  prescribed  by  the  lubrication 
theory  and  may  be  completely  solved  to  yield,  not 
only  the  film  pressure  distribution,  but  also  the 
normal  film  flux  at  all  points  of  the  film 
perimeter. 

While  the  film  continuity  condition  Is  Inher¬ 
ent  In  the  Interior  of  the  bearing  film  as  dic¬ 
tated  by  the  lubrication  theory,  the  transboundary 
continuity  condition  is  more  subtle.  One  may  bor¬ 
row  from  gas  dynamics  the  notion  of  the  Ranklne- 
Hugoniot  conditions  to  allow  for  abrupt  changes  of 
state  properties  while  preserving  continuity  of 
lubricant  transport.  As  Ranklne-Hugonlot  jumps 
are  approximations  of  the  diffused  transport  proc¬ 
esses  across  a  moving  shock  wave,  the  analogy  for 
a  moving  cavitation  boundary  is  the  change  of  flow 
structure  that  takes  place  within  a  distance  of 
the  scale  of  the  bearing  gap. 
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SOME  EXPERIMENTAL  FACTS 


Before  attempting  to  examine  the  equivalent 
Ranklne-Hugonlot  problem  at  the  Internal  bounaary 
of  an  incomplete  film,  It  Is  Instructive  to  recall 
some  representative  experimental  facts  of  Incom¬ 
plete  films.  Two  photographs  (Fig.  1)  which  first 
appeared  In  the  first  Leeds-Lyon  symposium  (29), 
reveal  that  the  notion  of  a  smooth  Internal  bound¬ 
ary  which  separates  a  filled  gap  from  a  void  re¬ 
gion  Is  at  best  an  approximation.  In  reality,  the 
void  region  has  a  textured  two-phase  flow  struc¬ 
ture.  The  border  between  the  filled  gap  and  a 
void  region  exhibits  Irregularities  associated 
with  the  two-phase  flow  structure.  To  circumvent 
such  complicated  details,  it  Is  convenient  tc  deal 
with  spatially  averaged  quantities  of  the  two-phase 
region  In  modeling. 

It  Is  clearly  revealed  that  a  portion  of  the 
void-dominated  region  Is  locally  filled  with  liquid 
bridging  both  surfaces.  Less  visually  evident,  but 
not  to  be  overlooked,  are  the  possibilities  that  a 
thin  liquid  layer  may  be  attached  to  either  but 
not  both  surfaces.  Thus,  a  complete  morphological 
description  of  the  void  dominated  region  must  In¬ 
clude  the  following  distinct  liquid  states: 

State  A  -  A  thin  layer  of  liquid  that  adheres 
to  the  sliding  (journal)  surface. 

State  B  -  A  liquid  streak  or  spot  that 
bridges  both  surfaces. 

State  C  -  A  thin  layer  of  liquid  that  clings 
to  the  stationary  surface. 

Excluded  from  present  consideration  are  minute  liq¬ 
uid  droplets  that  move  freely  between  the  two  sur¬ 
faces. 


Fig.  1.  -  Appearance  of  high  void  content  flow 
structures  (29).  (a)  Divergent  gap  between 
rotating  cylinder  near  stationary  plane,  (b)  Sep¬ 
aration  of  sphere  from  submerged  stationary  plane. 


EQUIVALENT  RANKINE-HUGONIOT 

The  equivalent  Rankine-Hugoniot  Is  concerned 
with  the  abrupt  changes  across  3  moving  Internal 
boundary.  It  deals  with  the  conservation  of  the 
liquid  phase  as  related  to  the  thin-fllm  flow  on 
one  side,  liquid-phase  transport  on  the  other 
side,  and  the  motion  of  the  boundary. 

Figure  2  shows  that  the  normal  unit  vector  of 
the  internal  boundary  at  E  Is  directed  into  the 
void  region  and  is  inclined  from  the  sliding  vec¬ 
tor  by  angle  a.  The  film  flux,  according  to  the 
lubrication  theory,  is 

4 ,  =»  -<12p)"'h3grad  p  +  1/2  Uh  [313 

In  the  void  dominated  region  liquid  content  is 
given  by  the  volumetric  fraction  or  the  fill  fac¬ 
tor  for  each  of  the  three  morphological  states; 
namely,  ©a,  ©b.  ©c-  The  velocity  of  the  interna! 
boundary  is  orthogonal  to  itself.  Therefore,  the 
"sweeping  flux  vector"  is 

Q  =  Vh { 1  -  <©A  +  9r  t  0C)}  [32] 

The  transvoid  liquid  transport  flux  is 

<|>  °  Uh  (9a  +  1/2  ©0)  C33] 

Conservation  of  the  liquid  phase  require; 


This  is  the  equivalent  Rankine-Hugoniot  condi¬ 
tion.  It  shows  that  the  sweeping  speed  of  the 
internal  boundary  depends  on  the  neighboring  film 
flux  and  the  adjacent  morphological  states.  A 
more  restrictive  form  of  this  equation  was 
originally  presented  by  Olsson  (8). 


IMPLICATIONS  OF  THE  EQUIVALENT  RANKINE-HUGONIOl 
CONDITION 


Equation  [363  applies  to  both  the  breakup 
boundary,  where  cos  a  >  0,  and  to  the  reformation 
boundary,  where  cos  a  <  0.  It  must  be  satisfied 
by  all  valid  solutions  of  incomplete  films. 

One  of  the  solution'  of  the  Olsson-Rankine- 
Hugoniot  condition  is  the  breakup  line  of  a  sta¬ 
tionary,  incomplete  film.  This  corresponds  to  the 
Stationary  breakup  point: 

grad  p  =  o  [37a] 

©A  +  1/2  0B  =  1/2  C 37b ] 

The  first  of  these  is  well  known,  but  the  second 
one  is  rarely  cited. 

Another  solution  of  the  stationary  incomplete 
film  refers  to  a  point  on  the  reformation  line  or 
Stationary  reformation  point: 


<4>-y-Q)»n  =  0 


[343 


-  (12p)“’h2{(cos  a)(8p/8x)  +  (sin  a)  (3p/3Z)} 


Substituting  Eqs.  [313  to  1333  inso  Eq.  [34], 
one  obtains 


+  U(co$  a) {1/2  -  <©A  +  1/2  ©B>}  =  0 


[383 


h2[(cos  a) Op/3x)  4  ssin  cJ<3p/3z)3 
+  U(cos  a)Li /?  -  <©A  +  i/2  ©B> 3 

-  V{1  -  (SA  *•  ©B  +  ©r>}  =  C  [353 

Rearranging,  one  finds 

vn  -  <eA  +  ©B  +  ec>3  » u  cos  «  [|  -  <©A  >■  |©B>] 

-  (12p)  ' h2(cos  a  ^  +  sin  a  [363 


Fig.  a.  -  Vector  diagram. 


This  relation  was  used  by  Jakobsson  and  Floberg 
(1)  to  find  the  shape  of  au  enclosed  cavity  of  a 
submerged  journal  bearing. 

For  s  moving  internal  boundary,  the  Olsson- 
Raiif’ne-liugonlct  condition  serves  as  a  link  be¬ 
tween  an  elliptic  field  of  the  fully  filled  domain 
and  a  hyperbolic  field  of  liquid  phase  transport. 

In  this  context  one  may  observe  that,  because  the 
normal  flux  is  a  unique  consequence  of  the  D1 rich- 
let  problem,  one  c?n  not  rule  out  the  possibili¬ 
ties  that  It  may  be  inwardly  directed,  and  that 
the  film  prassure  may  be  less  than  the  boundary 
value. 

UNRESOLVED  ISSUES 

In  order  to  use  the  O'sson-Rankine-Hugonlot 
conr  t'cr.  the  three  morphological  state  parame¬ 
ter  6*.  ©b>  anf1  eC  rau  t  i,e  9'vsn.  The  ©c 
probably  r  be  neglected  unless  there  are  special 
.h  .lica!  cdjses  for  a  bleed  affinity  of  the  lubri¬ 
cant  to  the  stationary  bear ^ 'k  surface.  Another 
p".vusU>lr.  hypothesis  is  to  hnpose  £q.  [37b]  at 
fch  b’-ea'rup  p..1nt,  irrespective  of  whether  the 
‘ntnrnal  nouiOa*  -  may  be  in  motion.  On?  is  thus 
lev.  with  the  remaining  question  of  how  liquid 
transport  beyond  the  meniscus  may  be  divided 
between  an  adhered  film  ©a  and  streamers  ©g. 

Sivt  .>se  that  they  are  gi'"T  initial  values  at 


-1 
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point  E;  then  subsequent  liquid-phase  transport 
due  to  sliding  can  be  stated  as 

<3/3t)[(0A  -  0B)h]  +  U(3/3x)[(0A  +  1/2  0B)h]  «  0 

[39] 

And,  if  one  is  willing  to  neglect  "cross  flowing" 
between  the  two  constituents,  the  streamer  trans¬ 
port  is  given  by 

(3/3t)(0Bh)  +  1/2  U(3/3x)(0Bh)  =  0  [39a] 

and  the  adhered  film  transport  Is  given  by 

(3/3t)(0Ah)  +  U  (3/3x)  <0Ah)  =  0  [39b] 


Clearly,  a  physically  permissible  initial  morpho¬ 
logical  condition  at  E  is  in  the  range 

0  3  <0A)E  S  1/2  [40] 

If  one  wishes  to  Improve  the  speculative  aspects 
of  these  ideas,  details  of  the  three-dimensional, 
viscous,  free  surface  flow  with  a  moving  surface 
and  a  moving  attachment  line  must  be  studied. 

The  implementation  of  the  Olsson-Rankine- 
Hugoniot  condition  in  a  numerical  computation 
scheme  is  not  trivial.  It  refers  to  the  film  flux 
at  the  boundary,  which  in  turn  shifts  with  time. 

It  is  quite  clear  that  the  computation  scheme  must 
be  capable  of  dealing  with  oblique  mesh  lines. 
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GEOMETRY  AND  MOTION  OF  THE  MENISCUS 


JEAN  FRENE,4  LUC  LAUNAY,4  AND  CODA  H.T.  PAN5 

Dynamic  loading  is  the  primary  concern  in  the 
design  for  journal  bearings  of  automotive  engines. 
Oil  supply  is  to  be  optimized  to  ensure  that  the 
minimum  gap  is  adequate  while  friction  is  to  be 
kept  as  low  as  possible.  Operation  with  an  incom¬ 
plete  film  is  typical  in  such  bearings. 

In  the  past  the  ir-film  model  is  commonly 
used  in  analytical  studies.  As  engineers  strive 
for  more  accurate  modeling  tools,  a  more  precise 
treatment  of  the  dynamics  of  the  incomplete  film 
is  in  order.  The  rr-fl lm  model ,  as  applied  to  the 
steady-state,  plain  journal  bearing,  is  equivalent 
to  the  Giimbel  approximation  in  rejecting  the  subam- 
Lient  portion  of  the  pressure  field.  In  conjunc¬ 
tion  with  the  short  bearing  approximation,  the 
ir-film  model  yields  simple,  closed-form  analytical 
results  for  both  static  (eccentricity)  and  dynamic 
(radial  velocity)  effects.  This  combination  of 
analytical  tools  has  been  the  backbone  for  the 
study  of  journal  bearings  in  automotive  engines. 

The  Reynolds-Swift-Stleber  cavitation  model  (RSS) 
has  been  regarded  to  be  more  precise.  Because  it 
requires  a  null  condition  of  the  pressure  gradient 
at  the  cavitation  boundary,  it  had  been  thought  to 
be  intrinsically  contradictory  to  the  short  bear¬ 
ing  approximation  and  had  been  used  primarily  in 
"FDM  type"  numerical  computations.  More  recently, 
a  modified  Reynolds-Swift-Stleber  model  was  used 
to  improve  the  short  bearing  approximation.  Cavi¬ 
ties  of  both  vented  and  enclosed  types  have  been 
successfully  modeled. 

The  Reynolds-Swift-Stieber  model,  which  is 
restricted  to  the  steady-state  condition  only,  is 
thus,  valid  for  either  a  static  load  or  a  rotating 
load.  As  it  was  applied  with  the  short  bearing 
approximation,  the  results  are  appropriate  for 
statid-load  considerations  of  the  journal  bearing 
or  for  high  rotating  loads  of  the  squeeze-film 
damper.  The  bearings  of  the  automotive  engine 
experience  oscillatory  loads  which  are  not  of  the 
pure  rotating  type,  and  the  static  load  component 
is  not  totally  negligible.  These  applications 
motivate  us  to  consider  the  more  complicated  prob¬ 
lem  of  a  truly  time-dependent  incomplete  film. 

The  short  bearing  approximation  was  retained  for 
this  work,  so  that  one  may  direct  full  attention 
to  issues  related  to  the  moving  cavitation 
boundary. 

ANALYTICAL  APPROACH 


The  time-dependent,  incomplete  film  problem 
requires  the  simultaneous  treatment  of  an  ellipti¬ 
cal  domain  and  a  hyperbolic  domain,  which  are 
joined  at  the  moving  internal  boundary  where  the 
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Olsson-Rankine-Hugoniot  condition  is  to  be  satis¬ 
fied.  A  centrally  fed  short  journal  bearing  will 
serve  as  the  specific  model  in  the  present  work. 
The  dimensionless  film  thickness  of  a  dynamically 
loaded  journal  bearing,  neglecting  angular  mis¬ 
alignment,  is 


H  =«  1  +  c  cos(9  -  4>)  [41  ] 


In  Eq.  [41]  both  e  and  4>,  may  be  time-dependent, 
employing  the  short  bearing  approximation  of 
Ocvirk,  the  Poiseuille  component  of  the  circumfer¬ 
ential  flux  is  neglected;  thus, 

~  1/2  wRCH  [42] 

¥z  -  -nzpHiooSop/az)  [43] 


The  continuity  condition,  with  allowance  for  mo¬ 
tion  of  the  journal ,  is 

(3<j/0/3e)  +  R(3»)»2/3z)  +  RC(3H/3t)  =  0  [44] 


Upon  making  use  of  Eqs.  [41]  and  [42],  one  finds 


(32p/3z2)  a  (12pC“2)H~3 


Sfcos  0  - 


- 


^jf-sin  0-<{) 


[45] 


Boundary  conditions  for  the  chosen  problem  are  of 
the  Dirichlet  type.  At  the  lubricant  supply  plane 

P(Z  =  0,0)  =  ps  [46] 

and  at  the  end  of  the  film 

p(Zg ,0)  =  0  [47] 

The  end  of  the  film  may  be  at  either  the  discharge 
end  of  the  bearing  or  at  the  internal  boundary  of 
the  incomplete  film;  i .e. , 

Zg  g  L/2  [48] 

The  elliptic  problem  can  thus  be  determined  by 
integrating  Eq.  [45]  with  the  constants  of  integra¬ 
tion  chosen  to  satisfy  Eqs.  [46]  and  [47],  to  yield 

p(z ,0)  a  ps(]  -2z/L) 

-  1/2  z(L/2  -  z)(12|i/C2)H"3[e  cos(0  -  <j>) 

-  (1/2  w  —  <j>)  sin(0  -  4>)]  [493 
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LIQUID  TRANSPORT  IN  THE  DISCHARGE  CAVITY 

In  the  discharge  cavity  liquid  phase  is 
assumed  to  be  in  the  state  of  an  adhered  film. 
The  liquid  content  is  given  by  the  adhered-fi 1m 
fill  factor  fA.  Transportation  of  the  adhered 
film  is  described  by  the  equation 


But  in  the  reformation  branch  fA  is  computed 
according  to  Eq.  C513,  so  that  Eq.  [52]  would  be 
used  in  its  fuli  form.  It  is  convenient  to  repre¬ 
sent  the  sweeping  motion  by 

(ft>e  =  «c  a  V 


<8/3t  +  u  0/30)  HA  =•  0  [50] 

where  HA  =  fAH.  The  time-stepped  representation 
of  this  solution  is 

HA(0,z;t)  =  Ha<0  -  <oAt ,  z;  t  -  At) 

=  1/2  H(0  —  cot  +  ut,  z;  t)  [51] 

for  t  >  t,  where  t  is  the  instant  the  adhered 
film  was  formed  at  the  breakup  point. 

MOTION  OF  THE  INTERNAL  BOUNDARY 

The  condition  of  01 sson-Rankine-Hugoniot 
gives  the  sweeping  speed  at  an  Internal  boundary 
to  be 

V  =  <1  ••  fA)_1  1/2  U(1  -  2fA)  cos  a 

-  <CH)2<12p)'"1(£2jpL||  +  sin  cc  ||)  [52] 

The  sweeping  speed  V  is  aimed  along  the  cavity, 
pointing  normal  to  the  periphery.  The  boundary  at 
any  instant  is  described  by 

ze<0,t)  3  L/2  [53] 

which  has  a  waist  at  0W  such  that 


=  (1  -  fA)-1[l/2  U  (1  -  2fA)  ctn  a 

-  (CH)2(12p)",(^|-^|g  +  sin  +  fa)]  [58] 

Then  updating  of  the  boundary  geometry  can  be 
executed  by  computing 

t 

zeC0,t)  =  ze(0,  t  -  At)  +  j  (3/3t)ze  dt  [59] 

t-At 

at  all  inner  points  of  the  internal  boundary  The 
end  points  move  along  z  =  L/2  and  can  be  calcula¬ 
ted  from 

(3/3t)x  =  sec  a  V 

e  z=L/2  z=L/2 

=  (1  -  fA)-1  1/2  U(1  -  2fA) 

-  *•»  «1)  z_L/2 

[60] 

Again,  fA  would  be  set  to  1/2  for  the  breakup 
point  but  would  be  computed  according  to  Eq.  [50] 
from  the  prior  history  of  the  reformation  point. 


<3/ae)z  L  =  0  [543 

e  °w 

The  inclination  of  the  periphery  normal  is 


Branches  on  either  side  of  the  waist  are 
joined  with  null  slope.  But  analytic  continuation 
from  one  to  the  other  across  the  waist  point  can 
not  be  assured  because  computation  of  Eq.  [51] 
would  follow  different  rules  for  the  two  branches. 
In  the  breakup  (0  <  9W),  the  Reynolds-Swift- 
Stieber  law  applies  so  that 

f A  =  1/2  [563 

and 

^breakup  "  -<C»)2<6p)-!(^2p-|  .  sin  a  ||)  [57] 


EQUATION  OF  MOTION 

The  journal  mass  may  experience  an  accelera¬ 
tion  as  the  bearing  film  force  may  not  cancel  the 
applied  force  completely.  This  acceleration  can 
be  integrated  in  time  to  determine  the  varying 
velocity  of  the  journal  center.  The  corresponding 
squeeze-film  terms  would  necessitate  recalculation 
of  the  fluid  film  forces  with  the  new  cavity 
boundary. 


PERIODICITY 

In  the  automotive  application  the  periodic 
state  of  the  constant  speed  operation  of  the 
engine  is  of  Interest.  The  character  of  the  above 
procedure  of  calculation  does  not  assure  this  peri¬ 
odicity  in  advance.  Therefore,  it  is  necessary  to 
perform  the  computation  for  more  than  a  fuil 
cycle,  until  the  results  demonstrate  complete 
attenuation  of  the  initial  transient  caused  by  the 
assumed  Initial  values  of  e  and  <j>. 
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CLASSIC  CAVITATION  MODELS  FOR  FINITE  ELEMENT  ANALYSIS 
J.F.  BOOKER8 


A  review  of  several  mathematical  models  com¬ 
monly  used  to  Incorporate  cavitation  In  the  finite- 
element  analysis  of  lubrication  is  presented.  The 
various  mathematical  models  are  completely  equiva¬ 
lent  In  the  limit,  but  they  differ  in  physical 
Interpretation  and  numerical  implementation.  The 
models  are  classic  In  the  sense  that  they  are  not 
flow  preserving.  For  clarity  we  review  the  devel¬ 
opment  of  finite-element  lubrication  analysis  both 
with  and  without  the  complication  of  cavitation. 

CAVITATION  EXCLUDED 

Transient  lubrication  analysis  Is  usually 
based  on  a  boundary-value  problem  in  which  the 
pressure  satisfies  the  Reynolds  differential  equa¬ 
tion  (30)  over  a  complete  film  region,  and  either 
pressure  value  or  normal  pressure  gradient,  respec¬ 
tively,  satisfy  "essential  ar.d  natural"  conditions 
everywhere  on  its  boundary,  whose  location  Is  fixed 
In  the  absence  of  cavitation.  In  an  equivalent 
variational  (extremum)  problem  (31),  a  pressure 
distribution  Is  sought  which  minimizes  a  func¬ 
tional  (Integral)  Involving  variants  of  the  Rey¬ 
nolds  operators  over  the  region  and  "essential" 
boundary  conditions  on  Its  boundary. 

Following  finite-element  discretization  of 
the  pressure  field,  the  functional  can  be  expressed 
In  terms  of  nodal  pressure  values.  Necessary  con¬ 
ditions  for  the  resulting  unconstrained  minimiza¬ 
tion  (quadratic  programming)  problem  give  a  linear 
algebraic  equation  set.  The  set  can  also  be  ob¬ 
tained  through  application  of  Galerkln's  method  to 
the  boundary  value  problem. 

Solution  of  the  linear  equation  set  com¬ 
monly  accomplished  by  either  direct  decomposition 
(e.g.,  Gauss,  Choleskl,  LUD  (lower  and  upper 
decomposition))  or  Indirect  Iteration  (e.g., 
Gauss-Seidel ,  and  SOR)  methods. 

CAVITATION  INCLUDED 

In  the  classical  view  of  transient  cavitation 
analysis,  the  complete  film  region  has  one  or  more 
free  boundaries  on  which  both  pressure  value  and 
normal  pressure  gradient  satisfy  so-called  Rey¬ 
nolds  boundary  conditions  (usually  null).  In  a 
second  view,  a  nonnegative  solution  (over  the  en¬ 
tire  film  region)  is  sought  for  the  variational 
problem.  The  resulting  constra  led  minimization 
(quadratic  programming)  problem  no  longer  admits  a 
simple  linear  algebraic  solution.  Third,  we  can 
choose  to  apply  a  mixture  law  view  of  two-phase 
(cavitating)  fluids  (15).  Applied  here,  a  com¬ 
plete  film  region  is  characterized  by  positive 
pressure  and  constant  (mixture)  density  value 
(unity  volume  fraction);  an  incomplete  (incipl- 
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ently  cavitating)  region  Is  characterized  by  null 
pressure  and  negative  density  time  rate  of  change. 
The  physical  Interpretation  is  obviously  that  of 
vaporization.  It  can  be  shown  that  the  three  view¬ 
points  are  equivalent  In  the  limit.  Any  one  is 
sufficient  In  principle  to  fix  the  (instantaneous) 
division  of  complete  and  incomplete  film  regions 
and/or  the  positions  of  their  (free)  boundaries. 

In  practice,  however,  all  three  viewpoints  Intro¬ 
duce  nonlinear  complexities  to  the  problem  and  a 
consequent  layer  of  Iteration  to  its  numerical 
solution. 

Focusing  on  boundary  location,  as  In  the 
first  view,  is  particularly  awkward  for  finite- 
element  analysis.  Since  locations  of  free  bounda¬ 
ries  are  unknown  a  priori,  they  must  be  selected  In 
iterative  trials  (with  attendant  nodal  discretiza¬ 
tion)  followed  by  solution  of  a  weli-posed  problem 
over  the  complete  film  region  (with  only  one  of 
the  two  Reynolds  boundary  conditions  fixed  and  the 
other  evaluated  for  iterative  closure). 

Focusing  on  local  behavior  Is  particularly 
convenient  for  finite-element  analysis,  since  It 
allows  a  fixed  nodal  discretization.  The  con¬ 
strained  quadratic  programming  problem  of  the  sec¬ 
ond  view  can  be  solved  In  a  number  of  ways,  all 
Iterative.  Physically  appealing  Is  Chrlstopher- 
son's  method,  In  which  the  linear  equation  set  Is 
solved  by  Gauss-Seidel  (or  SOR)  Iteration,  after 
which  any  computed  negative  pressures  are  set  to 
zero,  and  the  process  repeated  until  satisfactory 
convergence  Is  achieved.  (One  can  Imagine  an 
analogous  physical  process  in  which  quantities  of 
liquid  are  converted  to  vapor.)  The  scheme  has 
been  used  with  success  (32,33);  Its  main  drawback 
Is  Its  reliance  on  the  Iterative  convergence  of 
real  variables.  An  alternative  to  the  constrained 
quadratic  programming  problem  Is  the  linear  comple- 
mentarily  problem  (equivalent  In  the  limit).  Its 
numerical  application  (34-36)  appears  Identical  to 
that  of  the  mixture  law  (37),  in  which  one  seeks 
to  partition  nodes  as  cavitating  (or  not),  depend¬ 
ing  on  whether  particular  equality  and  inequality 
constraints  are  satisfied. 

A  metaphor  would  be  a  "porcupine  bearing"  In 
which  all  interior  nodes  communicate  via  feed/bleed 
check-valves  to  an  ambient  pressure  source.  Al¬ 
though  no  cavitation  at  all  would  occur  in  such  a 
preposterous  arrangement,  the  mathematical  problem 
is  Identical  to  those  above. 

Solution  of  all  such  problems  Is  a  matter  of 
integer  programming  (guesswork)  combined  with  di¬ 
rect  solution  of  linear  equation  sets  for  pres¬ 
sures  and  direct  evaluation  for  flow  terms.  In 
Murty's  method  (analogous  to  Christopherson's 
method)  nodes  at  which  constraint  violations  occur 
are  shifted  to  the  opposite  partition  for  the  next 
Iteration.  The  great  advantage  of  such  Integer 
iteration  is  that  convergence  is  abrupt  and  unam- 
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biguous;  with  a  bit  of  luck,  It  is  also  fast.  Fig¬ 
ure  3  shows  pressure  distributions  from  a  series 
of  such  steps,  starting  with  the  assumption  of  no 
cavitated  nodes  whatsoever.  It  should  be  noted 
that  such  a  blind  start  Is  unnecessary  In  practice; 
reasonable  first  Iterates  result  in  convergence  In 
one  or  two  steps.  (Corresponding  nodal  flow  quan¬ 
tities  do  not  provide  for  a  meaningful  plot  In 
this  previously  unpublished  study  by  LaBouff.) 

SUMMARY 

Several  mathematical  models  commonly  used  to 
incorporate  cavitation  in  the  finite-element  analy¬ 
sis  of  lubrication  and  suggested  both  physical  in¬ 
terpretation  and  numerical  implementation.  These 
classic  models,  which  are  equivalent  in  the  limit 
of  discretization,  describe  incipient  (quasis¬ 
tatic), nonevolving  (flow-preserving)  transient  cav¬ 
itation,  in  the  sense  that  vapor  only  begins  to 
form  before  the  liquid  is  replenished.  Though 
they  are  physically  unrealistic  in  a  really  most 
basic  way,  the  fact  that  they  fix  the  cavitation 
situation  quasistatical ly  in  instantaneous 
response  to  bearing  kinematics  means  that  the  de¬ 
grees  of  freedom  of  the  problem  are  hundreds  of 
times  fewer  than  if  nodal  density  (volume  frac¬ 
tion)  must  be  called  on  to  reflect  evolutionary 
history.  The  models  thus  lead  to  computationally 
efficient  algorithms  and  are,  therefore,  likely  to 
persist  for  a  long  time  to  come,  whatever  their  de¬ 
partures  from  realism. 
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MATHEMATICAL  MODELS  AND  RELATED  ALGORITHMS  IN  CAVITATION  PROBLEMS 


G.  BAYADA,7  M.  CHAMBAT,4  M.  EL  ALAOUP 

We  concern  ourselves  with  gaseous  cavitation 
operating  under  both  the  steady-state  and  transient 
conditions.  Applications  to  circumferentially  sup¬ 
plied  bearings  and  to  seals  are  given.  We  describe 
first  the  widely  used  cavitation  algorithms  and  the 
related  mathematical  problems  especially  in  view  of 
the  output  flow  computation.  We  propose  new  algo¬ 
rithms  using  the  finite  elements  variational  ap¬ 
proach  whose  reliability  has  been  theoretically 
proved.  For  unsteady  problems,  a  time  discretiza¬ 
tion  procedure  using  the  characteristic  lines 
allows  us  to  Introduce  a  set  of  stationary  like 
problems  for  which  equality  between  input  and  out¬ 
put  flows  is  preserved. 

STEADY  STATE  PROBLEMS 

Classical  Models 

The  Giimbel  (or  Half-Sommerfeld)  model  of 
cavitation  is  widely  used  because  it  is  easy  to 
Implement,  but  it  is  associated  to  a  false  bound¬ 
ary  condition  on  the  free  boundary  surface.  The 
situation  is  much  more  confusing  for  the  Swlft- 
Stieber  or  Reynolds  models.  The  validity  of  this 
last  model  has  been  supported  in  various  ways: 

a.  The  regularity  condition  is  required  to 
hold  at  the  free  rupture  boundary  (38): 

p  =  3p/3n  =  0  [61] 

b.  In  numerous  papers  condition  [61]  is  gener¬ 
alized  to  any  free  boundary  (39). 

c.  In  a  more  mathematical  formulation,  condi¬ 
tion  (61)  is  deduced  from  a  minimization  argument 
(40):  Through  all  positive  smooth  functions  satis¬ 
fying  the  external  boundary  conditions,  the  pres¬ 
sure  is  the  function  that  minimizes  the  energy 
integral . 

d.  The  pressure  is  the  solution  of  a  varia¬ 
tional  inequality  (41,42). 

e.  The  pressure  is  the  solution  of  a  comple¬ 
mentary  problem  (43). 

f.  The  pressure  is  obtained  by  the  Christoph- 
erson  procedure. 

There  are  two  types  of  unpleasant  problems 
that  appear  both  analytically  or  numerically: 
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1.  For  a  journal  bearing  with  an  axial  supply 
line  located  in  the  convergent  part  of  the  gap,  a 
strange  pressure  distribution  has  been  obtained 
numerically  with  nonconnected  active  areas  (44). 

2.  The  balance  between  the  input  flow  and  out¬ 
put  flow  is  not  always  satisfied,  especially  for 
seals  (45). 

Theoretical  studies  lead  to  the  following  conclu¬ 
sions  (46): 

Formulations  c,  d,  e,  f  are  equivalent.  They 
are  associated  with  a  well-posed  problem 
which  has  a  unique  solution  that  fulfills  the 
conditions  of  Eq.  [61].  This  solution  is 
such  that  the  cavitation  area  is  located 
only  in  the  divergent  parv  of  the  mechanism. 

A  New  Model 

Various  interpretations  of  the  model  are 
described  in  references  (1)  and  (14).  A  summary 
of  those  interpretations  is  given  here.  The  pres¬ 
sure  is  not  a  unique  parameter  governing  the  cavi¬ 
tation  phenomenon;  and  we  must  Introduce  another 
parameter  6  which  can  be  viewed  as  a  saturation 
function  or  a  fractional  film  content  (the  percent¬ 
age  of  oil  contained  In  an  elementary  volume)  or 
the  result  of  the  homogenization  of  the  oil  fin¬ 
gers.  In  any  case,  the  cavitated  area  is  not 
empty  as  It  contains  some  lubricant. 

The  flow  is  now  defined  by 

$  =  h3/(12p)Vp  -  6U  h/2  [62] 

assuming  (p  §  0  and  0  <  0  <  1),  with  the  con¬ 
straint  that  (0  »  I  If  p  >  0)  then  the  condition 
div  ($)  =0,  in  a  weak  variational  sense,  gives 
the  equations  for  the  problem. 

Numerical  Method  and  Applications 

Elrod  proposed  an  iterative  procedure  (15) 
well  adapted  for  finite  difference  discretization. 
We  propose  a  new  procedure  using  the  finite-element 
discretization,  whose  accuracy  and  convergence  have 
been  rigorously  proved  (46). 

It  is  possible  to  express  the  problem  in  terms  of 
one  single  unknown  X  which  is  a  solution  of  a 
nonlinear  equation 

X  =  B(X)  [63] 

This  equation  can  be  solved  using  the  fixed-point 
method  starting  from  zero.  r>nce  X  is  known,  it 
is  easy  to  recover  the  Initial  functions  p  and 
0. 

We  consider  a  bearing  with  a  circumferential 
supply  line.  We  compute  the  load,  the  attitude 
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angle,  and  the  input/output  flow  rates  for  various 
supply  pressures,  various  eccentricities,  and 
three  length  to  diameter  ratios.  The  first  compu¬ 
tation  is  related  to  the  Reynolds  model  solved 
using  the  Chrlstopherson  numerical  procedure,  and 
the  second  one  to  the  Elrod-Jacobson-Floberg  model 
using  the  new  procedure.  Both  computations  in¬ 
volved  finite  elements.  Results  demonstrate  lit¬ 
tle  discrepancy  for  the  load  between  the  two 
models,  but  gave  more  discrepancy  for  the  eccen¬ 
tricity  angle.  The  equality  between  output  flow 
and  Input  flow  is  very  well  preserved  for  the  sec¬ 
ond  model.  However,  it  is  not  true  for  the  first 
model  for  which  the  discrepancy  Increases  with 
eccentricity  and  length  to  diameter  ratio. 

TRANSIENT  PROBLEMS 


Classical  Models 


The  same  analysis  can  be  made  for  steady- 
state  problems  from  both  the  theoretical  and 
numerical  aspects.  The  Chrlstopherson  procedure 
is  commonly  used  with  an  explicit  or  Implicit  time 
descretization.  Bad  results  are  usually  reported 
for  the  balance  of  the  flow,  especially  for 
seals.  These  will  be  Illustrated  as  a  prime 
example. 


3(6h)  a/hi.  8fi\  1  3(ehU) 
3t  “  3h\J2p  dxj  +  2  3x 


[64] 


This  equation  Includes  all  the  above  equations, 
especially  the  condition  on  the  free  boundary, 
which  appears  exactly  as  the  classical  Ranklne 
Hugonlot  condition  for  shocks  In  hyperbolic  equa¬ 
tions  and  is  equivalent  to  those  proposed  by  Elrod 
(14).  Supplementary  data  involve  (1)  boundary 
conditions  on  the  pressure  and  (2)  initial  data 
for  0. 

Numerical  Procedure  (47) 


The  starting  point  Is  to  substitute  the  total 
derivative  0(6h)/Dt  for  l/2C3(0hu)/3x]  +  (36h/3t) 
with  the  meaning  that  the  derivative  in  t  is  taken 
along  the  direction  U  so  that  the  equation 
becomes 


D(6h)  Lf^So)  n 
Dt  "  3x^l2}i  3 xj 


0  =  1  if  p  >  0  C65] 


and  then  discretize  the  first  term  of  the  equation 
in  time  along  the  characteristic  line. 

CONCLUSIONS 


Mew  Model 

Me  describe  It  first  In  the  particular  case 
of  an  infinitely  long  mechanism.  However,  the  gen¬ 
eralization  to  real  two-dimensional  geometry  is 
straightforward. 


An  application  to  a  seal  was  made  in  (47). 
Rectangular  finite  elements  were  used. 

1.  The  discrepancy  between  input  mass  flow 
and  output  mass  flow  never  exceeds  4  percent  In 
each  period. 


As  In  the  steady-state  problem,  a  unique 
equation  in  a  weak  or  distributional  sense  Is  well 
suited  for  numerical  procedure: 


2.  Periodicity  of  the  load  and  of  the  flows 
Is  obtained  In  four  periods.  The  cavity  follows 
the  rotor  movement  when  X£  =  2xj. 
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NON-NEWTONIAN  EFFECTS  IN  FILM  CAVITATION 


A.Z.  SZERI 1 0 ,  K.R.  RAJAGOPAL10,  AND  ZHANG10 

There  exist  in  engineering  usage  a  vast 
number  of  fluids,  both  naturally  occurring  and  man 
made,  the  behavior  of  which  cannot  adequately  be 
described  by  classical,  linearly  viscous  theory. 
Though  such  fluids  also  abound  in  tribological 
practice,  we  have  accumulated  only  scant  knowledge 
of  the  effects  non-Newtonian  behavior  might  have 
on  tribological  performance.  The  principal  prob¬ 
lem  we  are  faced  with  In  the  area  of  non-Newtonian 
fluid  mechanics  Is  our  Inability  to  successfully 
model  the  stress-deformation  relationship  of  a  par¬ 
ticular  material.  But,  In  part,  we  can  also  blame 
the  lack  of  reliable  experimental  data  and,  thus, 
the  absence  of  guidance  to  theoretical  research. 

Viscoelastic  fluids  differ  from  Newtonian 
fluids  In  one  main  aspect,  viz.,  viscoelastic 
fluids  possess  memory.  Having  memory,  for  the 
fluid,  means  that  the  stress  presently  experienced 
by  a  fluid  particle  Is  no  longer  determined  by  ex¬ 
isting  conditions  alone  but  depends  on  deformation 
history  of  the  material  In  the  neighborhood  of 
that  particle.  The  principal  practical  manifesta¬ 
tions  of  non-Newtonian  behavior  are  normal  stress 
effects  and  shear  thinning.  We  have  recently 
noted  another  strong  effect  on  the  velocity  pro¬ 
file  In  a  thrust  bearing-like  geometry  (48) 

(Fig.  4).  Velocity  profiles  were  measured  using 
laser  Doppler  .tnemometry,  In  a  flow  between  rota¬ 
ting  parallel  disks.  Figure  4  displays  the  tangen¬ 
tial  velocity  field  at  a  fixed  radial  position, 
with  the  top  disk  stationary  and  the  bottom  disk 
rotated  at  constant  angular  velocity.  The  profile 
In  (a)  was  obtained  In  water,  and  that  In  part  (b) 
was  measured  In  a  viscoelastic  fluid  (0.023  per¬ 
cent  polyacrylamide  solution).  It  Is  easy  to  dis¬ 
cern  the  qualitative  difference  between  the  two 
profiles:  Shear  thinning  occurs  at  a  certain  axi¬ 
al  position,  the  separation  distance  of  which  from 
the  bottom  plate  depends  on  various  parameters,  re¬ 
sulting  in  large  velocity  gradients.  The  cause 
for  this  localized  action  is  not  known  to  us. 

We  have  also  noted  recently  (49-51)  that  jour¬ 
nal  bearings  lubricated  with  non-Newtonian  lubri¬ 
cants  could  be  less  susceptible  to  thermal  effects 
when  compared  with  bearings  lubricated  with  Newto¬ 
nian  fluids. 

In  this  note  we  wish  to  further  explore  lubri¬ 
cation  with  viscoelastic  fluids  and  focus  our  at¬ 
tention  on  cavitation,  first  wishing  to  place  the 
constitutive  relationship  of  the  class  of  fluids 
we  employ,  In  proper  perspective.  To  this  end  we 
define  the  motion  of  a  fluid  particle  as  follows. 

If  fluid  particle  X  at  time  t  occupies  the  spa¬ 
tial  position  x,  we  write 

x  =  x(X,  t)  t  >  0  C66 ] 
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As  t  varies,  so  will  vary  the  position  x  of 
the  particle  X.  We  Insist  that  this  relationship 
(Eq.  [66])  be  invertible;  l.e.,  that  two  particles 
do  not  at  the  same  time  occupy  the  same  spatial 
position,  and  write  for  the  inverse  of  Eq.  [66] 

X  =  x_1(x,t)  [67] 

This  Is  just  mathematical  shorthand  for  the  state¬ 
ment  that  spatial  position  x  at  time  t  nas  the 
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fluid  particle  X  as  its  temporary  resident. 

Consider  now  another  time  Instant  x  <  t,  at 
which  time  our  fluid  particle  X  was  resident  at 
the  position  ?.  Hence,  we  have 

%  =  x(X,x)  [68] 

But  we  can  substitute  for  X  from  Eq.  (67)  into 
Eq.  [68]  and  have 

5  =  xCx"1<x,t),x]  =  xt^x ,x)  [69] 

The  relationship  §  =  xt(x,t)  states  that  the  fluid 
particle  that  is  in  position  x  at  the  present 
time  t,  was  in  position  %  at  some  former  time 
instant  we  call  x.  If  we  now  let  x  §  0  vary 
while  calling  t  the  present  time  and  keeping  it 
fixed,  we  are  looking  backwards  in  time  and  trac¬ 
ing  the  history  of  motion  of  the  fluid  particle  X. 

The  vector  valued  function  xt  occurring  in 
Eq.  [683  is  called  the  relative  deformation  func¬ 
tion.  If  the  deformation  xt  described  is  smooth 
and  the  material  has  no  microstructure  (and  this 
applies  to  most  deformations  and  most  materials  in 
engineering  practice),  we  can  define  the  gradient 
of  xt>  called  F(s),  through  (52). 

dS  =  F(s)  dx;  F(s)  =  9S/9x;  s  =  t  -  x  [70] 

Equation  [70]  shows  that  the  gradient  of  xt>  in 
our  notation  F(s),  Is  indicative  of  how  an  Infini¬ 
tesimal  line  element  d?  at  the  instant  of  time 
x  had  to  be  stretched  and  rotated  so  as  to  become 
the  line  element  dx  of  the  present  time  t.  As 
we  let  s  >  0  increase,  the  equality  here  signi¬ 
fying  present  time,  we  have  the  history  of  deforma¬ 
tion  characterized  by  F(s). 

The  deformation  gradient  history  F(s)  is 
used  in  the  definition  of  a  class  of  fluids  called 
"simple  fluids".  Simple  fluids  (and  the  adjective 
"simple"  refers  to  the  lack  of  microstructure  and 
does  not  imply  that  the  fluid  is  elementary  In  any 
sense)  are  those  in  which  the  stress  T  at  a 
point  depends  on  the  history  of  deformation  of  a 
neighborhood  of  that  point,  or 

s=® 

T  =  f  (F(s))  [713 

s=0 

There  is  a  fairly  good  agreement  among  practition¬ 
ers  of  non-Newtonian  fluid  mechanics  concerning 
the  validity  of  Eq.  [71]  and  the  assumptions  lead¬ 
ing  up  to  it.  But,  there  is  disagreement  on  how 
to  proceed  further,  for  Eq.  [70]  is  hardly  a 
practical ly  useful  statement. 

There  are  many  ways  to  approximate  the  func¬ 
tional  f  in  Eq.[71].  Integral  representations 
are  well  received  by  the  rheology  community  in  gen¬ 
eral  but  have  to  be  as  yet  handled  numerically  in 
a  practical  application.  Differential  representa¬ 
tions,  on  the  other  hand,  yield  to  numerical  treat 
ment  but  allow  only  for  infinitesimal  memory. 


It  has  been  shown  that  under  certain  condi¬ 
tions  Eq.  [71]  has  an  expansion  in  terms  of  cer¬ 
tain  differential  quantities.  When  this  expansion 
is  terminated  at  a  low  order  and  when  account  is 
made  of  the  particular  geometry  common  to  most 
lubricant  films,  we  essentially  obtain  that  the 
stress  has  the  following  structure  (cf.  (50)) 

T  =  -  pT  +  [p  +  |3(tr  A,2)]  A,  +  oA2  [72] 

Here,  p  is  the  molecular  viscosity,  a  and  |3 
are  non-Newtonian  parameters,  and  Aj  and  A2  are 
the  first  and  second  Rivlin-Ericksen  tensors, 
respectively,  defined  by 

♦  +  T 

Aj  =  (grad  v)  +  (grad  v) 

[73] 

~  ~  ~  -»  ->  t~ 

A2  =  dAj/dt  +  A^(grad  v)  +  (grad  v)  A^ 

The  material  parameter  a  promotes  normal  stress 
effects  while  a  (3  <  0  portrays  shear  thinning 
and  a  |3  >  0  will  make  the  fluid  shear  thicken. 
This  is  easily  seen  when  Eq.  [72]  is  written  for 
simple  shear  flow,  for  then 

xxy  =  [p  +  |3(du/dy)2](du/dy)  [741 

Models  of  this  type  (Eq.  [72])  h  e  been  used  by 
several  researchers  (53-55).  We  have  applied  it 
to  journal  bearings  (50,51)  and  demonstrated  that 
even  a  slight  departure  from  Newtonion  behavior, 
as  measured  by  the  numerical  values  of  a  and  (3, 
renders  bearing  performance  relatively  Insensitive 
to  changes  in  lubricant  temperature. 

In  a  continuation  of  previous  effort,  we 
looked  at  non-Newtonian  effects  on  cavitation  in  a 
journal  bearing.  Figure  5  shows  the  slow,  gradual 
change  in  the  pressure  profile  as  the  material 
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Fig.  5.  -  Pressure  profiles  in  shear  thickening 
fluid  (y  =  -0.1). 
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Fig.  6.  -  Pressure  profiles  In  shear  thinning 
fluid  (Y  =  -0.1). 


Fig.  7.  -  Foldlines  (S)  on  solution  manifolds  (M). 

parameter  3  Is  made  progressively  smaller  while 
a  is  kept  constant.  To  be  precise,  It  Is  not  a 
and  p  that  are  used  as  the  parameters  of  the  prob¬ 
lem  but  their  nondlmenslonal  counterparts  y  and 
8,  defined  by 

y  =  —  5  =  ^  C75] 

pR  p(T 

Here,  U  Is  the  surface  velocity  of  a  shaft  of 
radius  R,  and  C  is  the  radial  clearance  In  the 
bearing. 

Figure  5  shows  pressure  profiles  for  shear 
thickening  fluid,  5  >  0,  at  y  =  -0.1.  As  the 
value  of  S  decreases  from  S  =  0.004,  the  maxi¬ 
mum  pressure,  gradually  decreases,  otherwise  the 
profile  seems  to  change  little. 

Figure  6  shows  pressure  profiles  again  at 
y  =  -0.1  but,  now,  for  a  shear  thinning  fluid. 
Starting  at  8  =  -0.004,  the  pressure  curve  shows 
no  cavitation.  Increasing  the  value  Oi  o  siowiy 
Increases  the  maximum  pressure,  but,  otherwise, 
the  shape  remains  Invariant,  still  showing  no  cavi¬ 
tation.  This  phenomenon  persists  up  to  the  value 
5  =  0.0025,  a  slightly  shear  thickening  fluid. 

But  at  a  somewhat  larger  value  of  8  =  0.005, 

Fig.  5  already  shows  cavitation. 
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Fig.  8.  -  Pressure  profiles  at  y  =  -0.1. 
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Fig.  9.  -  Pressure  maxima  at  y  a  -0.1. 


Our  Investigations  are  Incomplete  at  this 
stage,  and  conclusions  are  not  well  founded.  But 
they  suggest  the  presence  of  a  fold  line  on  the 
solution  manifold  at  0.0025  <  8  <  0.005  when 
y  =  -0.1.  Here,  the  two  solutions,  one  with  cavi¬ 
tation  and  one  without,  join.  Of  course  It  Is  pos¬ 
sible,  or  even  likely,  that  the  noncavltating  solu¬ 
tion,  If  It  at  all  exists,  represents  an  unstable 
branch  and  is  therefore  not  reproducible  In  the 
laboratory.  All  we  suggest  at  this  juncture  Is, 
that  there  might  be  an  other  than  the  classical, 
cavltatlng  solution  to  the  problem.  Our  next  task 
Is  to  trace  the  fold  line  on  the  solution  manifold 
using  continuation  methods.  Figure  7  Illustrates 
some  possible  fold  lines  of  an  equilibrium  surface. 

Figure  8  shows  two  pressure  profiles,  one 
obtained  with  3  =  0.0025  exhibiting  noncavita¬ 
tion;  while  the  other,  calculated  with  3  »  0.005, 
exhibiting  a  classical,  cavltatlng  profile. 

Except  for  the  presence  or  absence  of  cavitation, 
the  profiles  are  virtually  identical.  Figure  9 
displays  the  pressure  maxima  at  y  =  -0.1  as  the 
parameter  3  Is  varied. 
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It  is  well  known  that  nonlinear  problems 
often  yield  a  rich  multiplicity  of  solutions,  sup¬ 
ported  on  identical  steady-state  boundary  condi¬ 
tions.  As  an  example  we  may  quote  the  work  of 
Bejamin  and  Mull  in  (56),  who  found  21  solutions  in 
the  flow  between  short,  concentric  rotating  cylin¬ 
ders.  Their  work  was  experimental  and  thus  shows 


much  more  than  idle  mathematical  curiosity.  This 
of  course,  constitutes  no  proof  that  what  we  have 
here  in  the  non-Newtonian  journal  bearing  is  also 
multiplicity  of  solutions;  but  it  perhaps  suggests 
that  our  findings  should  not  be  discarded  out  of 
hand. 
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CAVITATION  MODELING  FOR  ANALYSIS  OF  FLUID  FILM  JOURNAL  BEARINGS 
L.E.  BARRETT"  AND  L.A.  BRANAGAN" 


A  reasonable  description  of  the  cavitation 
phenomenon  in  thin  fluid  films  is  critical  to  the 
accurate  prediction  of  the  pressure  distribution, 
and  hence  forces,  in  fluid-film  journal  bearings. 
Cavitation  occurs  when  the  fluid  film  is  subjected 
to  conditions  of  tension  which  tend  to  produce  neg¬ 
ative  pressures  within  the  film.  Usually,  under 
these  circumstances,  the  pressure  falls  to  essen¬ 
tially  the  vapor  pressure  of  the  lubricant,  but  no 
lower.  Gaseous  components  of  the  lubricant  come 
out  of  solution  and,  together  with  vaporized  lubri¬ 
cant,  form  a  two-phase  flow  field  of  bubbles  and 
streamers.  Cavitation  usually  occurs  in  regions 
of  diverging  film  thickness  within  the  bearing  and 
depends  on  the  supply  pressure  of  the  lubricant  to 
the  bearing.  Hence,  cavitation  regions  can  be 
drastically  altered  by  changes  in  bearing  geometry 
and  supply  pressure,  and  for  a  given  bearing  they 
are  dependent  on  operating  parameters  such  as 
shaft  speed  and  bearing  load.  Figure  10  shows  a 
comparison  by  Tonnesen  of  predicted  and  measured 
pressures  obtained  in  a  two-axial -groove  journal 
bearing.  This  figure  clearly  shows  the  cavltated 
region  of  low  pressure  that  occurs  in  the  top  pad 
of  the  bearing  In  which  a  diverging  film  thickness 
occurs.  It  should  be  noted  that  cavitation  within 
bearing  fluid  films  is  not  considered  to  be  inher¬ 
ently  bad.  Only  rarely  has  cavitation  been  cited 
as  the  primary  cause  of  bearing  damage  or  failure 
due  to  surface  pitting,  and  cavitation  control  is 
rarely  a  goal  of  bearing  design  except  in  the  case 
of  squeeze-film  bearings. 

The  fluid  mechanics  and  chemistry  of  cavita¬ 
tion  Is  very  complex.  Fortunately,  for  most 
practical  design  considerations,  a  detailed 
description  of  the  pnenomenon  is  not  required. 

What  Is  required  is  an  accurate  prediction  of  the 
cavitation  region  boundaries  and  the  pressure,  tem¬ 
perature,  and  shear  heat  generation  within  the 
region.  This  is  usually  accomplished  by  applying 
reasonable  auxiliary  conditions  on  the  fluid  film 
pressures  and  temperatures  obtained  from  a  coupled 
iterative  solution  of  the  governing  equations,  the 
Reynolds  equation  for  pressure,  and  the  energy 
equation  for  temperature.  The  solutions  are  depen¬ 
dent  on  boundary  conditions  of  supply  pressure  and 
temperatures  and  may  include  the  heat  conduction 
effects  in  the  bearing  and  shaft.  The  simplest  so¬ 
lution  is  to  solve  the  governing  equations  indepen¬ 
dent  of  the  existence  of  cavitation  and  to  set  any 
predicted  negative  pressures  to  the  cavitation 
pressure  using  a  forward  solution  for  the  pressure 
starting  at  the  leading  edge  of  the  film.  It  is 
easily  shown  that  this  approach  violates  mass  flow 
continuity,  but  it  does  provide  a  reasonable  ap¬ 
proach  to  defining  the  cavitation  region  for  many 
bearing  applications,  particularly  if  certain  sim- 


"Oepartment  of  Mechanical  and  Aerospace  Engi¬ 
neering,  School  of  Engineering  and  Applied  Sci¬ 
ence,  University  of  Virginia,  Charlottesville,  VA. 


plifying  approximations  to  Reynolds  equation  are 
used,  such  as  the  short  bearing  approximation.  It 
is  called  the  half-Sommerfeld  pressure  boundary 
condition. 

For  finite-length  solutions,  this  approach 
has  usually  been  used  to  avoid  an  iterative  solu¬ 
tion  for  the  trailing  boundary  of  cavitated  region 
where  the  pressure  is  the  cavitation  pressure  and 
the  pressure  gradient  is  zero,  thus  ensuring  mass 
flow  conservation.  This  is  referred  to  as  the 
Reynolds  boundary  condition.  Actually,  the  trail¬ 
ing  cavitation  boundary  can  be  easily  found  with¬ 
out  iteration  by  back  solving  for  the  pressure, 
starting  from  the  trailing  edge  of  the  fluid  film 
and  substituting  the  cavitation  pressure  whenever 
the  solution  predicts  a  subcavitation  pressure. 

If  coupled  pressure-temperature  solutions  are 
being  sought,  the  problem  is  already  highly  Itera¬ 
tive,  so  there  is  no  advantage  to  be  gained  by 
using  the  half-Sommerfeld  cavitation  boundary  con¬ 
dition,  even  if  a  forward  substitution  scheme  is 
employed. 

Once  the  cavitation  boundaries  are  known,  the 
heat  generation  within  the  cavltated  region  must 
be  calculated.  The  heat  generated  has  an  indirect 
effect  on  the  cavitation  boundaries  and  no  effect 
on  the  pressure  within  the  cavltated  region.  The 
effect  on  the  boundaries  comes  about  because  the 
heat  generation  alters  the  temperature  within  the 
fluid  film  in  the  cavltated  region.  Thus,  the  tem¬ 
perature  of  the  lubricant  exiting  the  trailing 
edge  of  the  film  is  altered.  In  turn,  it  affects 
the  temperature  of  the  film  at  the  leading  edge  of 
the  following  pad. 

The  usual  theoretical  approach  is  to  use  what 
is  called  hot  oil  carry  over  calculations.  In 
this  method  the  lubricant  exiting  the  trailing 


Fig.  10.  -  Predicted  and  measured  pressures  in  a 
two-axial -groove  bearing. 
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Fig.  11.  -  Predicted  and  measured  pressures  in  a 
two-axial-groove  bearing.  Experiemnt 
conditions:  Tonnesen  bearing;  3500  rpm;  5600  N; 
Ts  =  50  °C.  Present  work  used  2°  elements  and  30 
harmonics . 


TOP  VIEW 
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SIDE  VIEW 


Fig.  12.  -  Model  for  two-phase  flow  in  cavitated 
lubricant  films. 

edge  of  a  pad  is  assumed  to  cling  to  the  shaft  and 
to  be  carried  into  the  Inlet  of  the  following 
pad.  If  this  volume  of  lubricant  is  insufficient 
to  fully  supply  the  following  pad,  the  remaining 
volume  Is  made  up  of  lubricant  supplied  through 
the  inlet  supply  groove.  The  temperature  of  the 
inlet  lubricant  is  determined  by  the  ratios  of  the 
carry  over  and  the  supply  groove  volumes  at  their 
respective  temperatures. 

A  detailed  theoretical  model  of  the  two-phase 
flow  in  the  cavitated  region  is  usually  not  con¬ 
ducted.  Some  assumption  of  the  flow  field  is 
therefore  made,  and  the  shear  heating  of  the  lubri¬ 
cant  1r.  this  region  is  determined  by  supplemental 
calculations.  The  simplest  model  is  to  assume 
that  the  lubricant  film  extends  from  the  bearing 
surface  to  the  pad  surface.  The  axial  film  width 
decreases  in  the  diverging  circumferential  direc¬ 
tion  since  there  is  a  constant  volume  of  lubricant 
and  because  there  is  no  pressure  gradient  deve¬ 


Fig.  13.  -  Predicted  and  measured  temperatures  in 
a  two  axial  groove  bearing  with  two-phase  flow 
model  used  in  cavitated  region.  Effect  of 
cross-film  bearing;  Tonnesen  bearing;  3500  rpm; 
Babbitt  temperature  profile;  1289.5  1 bf  load;  w/ 
thermal  cav. 

loped  in  the  cavitated  region.  This  model  treats 
the  cavitated  region  as  one  with  a  central  streamer 
of  lubricant  flowing  along  the  bearing  centerline. 
The  heat  generation  is  easily  calculated  since  the 
velocity  gradient  across  the  film  is  essentially 
constant,  at  least  for  laminar  flow,  being  approxi¬ 
mately  Rw/2,  v/here  R  is  the  shaft  radius  and 
u  is  the  shaft  angular  velocity. 

Calculations  performed  using  this  model  are 
shown  in  Fig.  11  for  a  two-axial -groove  bearing. 
This  figure  shows  the  circumferential  temperature 
variation  of  the  bearing  surface  in  both  the 
loaded  pad  and  cavitated  pad.  Measured  tempera¬ 
tures  made  by  Tonr.esen  are  shown.  Although  heat 
conduction  into  the  shaft  and  bearing  surface  is 
Included  in  the  theoretical  calculations,  the 
analysis  predicts  significantly  higher  tempera¬ 
tures  in  the  cavitated  region  than  those  meas¬ 
ured.  This  suggests  that  the  heat  generation 
model  is  predicting  larger  shear  heating  than 
actually  occurs. 

Figure  12  shows  a  two-phase  lubricant-gas 
model  for  heat  transfer  in  the  cavitated  region  in 
which  the  lubricant  adherence  to  the  shaft  surface 
can  be  assumed.  Further,  it  is  assumed  that  the 
height  of  the  film  diminishes  along  the  diverging 
portion  of  the  gap  rather  than  have  the  film 
decrease  in  width.  This  model  has  been  developed 
by  observations  of  flow  fields  in  transparent  bear¬ 
ings,  which  suggest  that  this  type  of  flow  pattern 
exists.  Because  the  viscosity  of  the  gaseous  por¬ 
tion  of  the  film  is  considerably  lower  than  that 
of  the  lubricant,  most  of  the  shearing  occurs  in 
the  gaseous  region.  Thus,  this  model  generates 
less  shear  heating  in  the  cavitated  region.  The 
heat  conduction  to  the  shaft  and  bearing  surfaces 
is  also  modified  because  of  the  difference  in 
thermal  properties  of  the  two  fields. 

Figure  13  shows  the  results  of  calculations 
performed  on  the  same  two-axial-groove  bearing  as 
presented  in  Fig.  11.  The  results  show  better  cor¬ 
relation  with  the  measured  temperatures  in  the  cav¬ 
itated  region,  suggesting  that  the  two-phase  model 
used  is  reasonably  good. 
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A  MODEL  FOR  CAVITATION  IN  DYNAMIC  BEARINGS 


JOHN  A.  TICHY1 2  AND  C.P.  KU'* 

In  fluid  film  lubrication  problems,  pressure 
often  seeks  a  level  lower  than  lubricant  vapor 
pressure  in  the  divergent  region  of  the  bearing, 

In  which  case  gas  or  vapor  forms  a  cavity.  For 
the  steadily  loaded  bearing,  the  cavity  Is  usually 
formed  due  to  ventilation  from  the  surrounding 
atmosphere  or  the  emission  of  dissolved  gases  from 
the  lubricant.  For  the  dynamically  loaded  bear¬ 
ing,  the  cavity  Is  composed  mostly  of  lubricant 
vapor. 

For  the  heavily  and  moderately  loaded  bear¬ 
ing,  the  Jakobsson-Floberg-Olsson  cavitation  model 
seems  to  be  correct  to  predict  the  extent  of  the 
cavitation  region  and  to  specify  proper  pressure 
boundary  conditions  on  the  Interface.  However, 
this  model  Is  difficult  to  apply  due  to  the  nonli¬ 
nearity  of  the  reformation  boundary  conditions. 

It  requires  extreme  numerical  effort  In  guessing 
the  shape  of  the  cavitation  region  and  checking 
the  flow  balance  (57). 


In  the  proposed  model  we  assume  the  shape  of 
the  pressure  axial  profile,  and  the  Floberg  bound¬ 
ary  conditions  are  maintained.  After  Integrating 
the  Reynolds  equation  across  the  axial  direction, 
a  nonlinear  ordinary  differential  equation  for  the 
distance  of  the  cavitation  boundary  from  the  bear¬ 
ing  centerline  as  a  function  of  the  circumferen¬ 
tial  coordinate  Is  obtained.  In  the  present  study 
we  consider  the  circular -centered  orbit  squeeze- 
film  damper  bearing  with  two  circumferential  oil 
supply  grooves  at  the  bearing  ends  (Fig.  14).  The 
lubricant  properties  and  temperature  are  assumed 
to  be  constant,  and  the  fluid  Is  laminar  and  iner¬ 
tialess.  The  unwrapped  dimensionless  bearing  sur¬ 
face  and  pressure  boundary  conditions  are  shown  in 
Fig.  15,  where  Z  =  0  Is  the  centerline  anu  W(X) 
is  the  width  of  the  cavity. 

ANALYSIS 


The  dimensionless  Reynolds  equation  and  the 
Floberg  boundary  conditions  are 
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where  C  Is  shown  In  Fig.  15  X  =  L'/Ro,  and  X], 
*2.  xm.  H(X)  are  unknowns.  The  P(X,Z)  Is  a  peri¬ 
odic  function  In  the  X  direction  and  global  mass 
conservation  requires 
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[78] 


After  Integrating  the  Reynolds  equation 
across  the  Z  direction  and  applying  boundary  con¬ 
ditions  at  Z  =  0  and  Z  »  1  or  H(X),  an  Integral 
Reynolds  equation  In  regions  I,  II,  and  III  Is 
obtained: 

Region  I: 


,,dH  lln 
12dX  +  ~x9zl 


[79] 


Fig.  15.  -  Bearing  surface. 
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Region  II: 

Jxf"3  f)  ’  ,2"  -  <  -  %.  (80> 

Region  III: 

5x("3  f)  -  '2<'  -  <  *  ^ZL  *  ,2[H<«>  -  HII 

.31] 

where  the  Integral  pressure  P(X)  is  defined  as 
the  Integral  of  P(X,Z)  across  the  Z  direction. 

Equations  [82]  and  [83]  are  pressure  profile 
assumptions  In  the  noncavitation  and  rupture  re¬ 
gions,  respectively.  These  pressure  profiles  sat¬ 
isfy  all  boundary  conditions  at  Z  =  0,1  and  Z  = 
W. 


PN<X,Z)  =  P$Z2  +  PC(X)<1  -  Z2>  [82] 


[83] 
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Fig.  16.  -  Cavity  shape. 
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|£<X,W)  +  ||(X,W)|£  =  0  [84] 

Calculating  the  Integral  pressure  of  Eqs.  [82]  and 
[83]  and  applying  the  chain  rule  (Eq.  [84]),  we 
substitute  these  equations  Inco  Eqs.  [79]  to  [81] 
and  apply  the  boundary  conditions  at  X],  Xm,  and 
X2  on  the  Interface  to  obtain  the  following  ordi¬ 
nary  differential  equations: 

Region  I: 


Fig.  17.  -  Centerline  pressure. 
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Region  III: 


P  =  |p$(l  -  W)  +  2X2[H(?)  -  H](  1  -  W)2  ^3  ^  [89] 
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Region  II: 


nVAi>  =  0 


LSI  3 


d2W  3  dH  dH  6 _ ]_ 

dx2  +  H  dX  dX  "  x2  1  -  W 


36(1  -  H)  dH  n 

+  <5  jy  a  0 

PSH3  dX 


W(X,)  =  ^(Xm)  =  0,  ^(X,)  =  DWI  [88] 


In  Eqs.  C85]  to  [91],  the  Independent  variable  Is 
X;  the  dependent  variables  In  the  three  regions 
are  Pc,  W,  and  P;  and  H  is  a  known  function. 
Known  parameters  are  Pg  and  X.  Unknown  parameters 
to  be  determined  from  the  solution  procedure  are 
the  locations  X],  Xm,  X2  (see  Fig.  15)  and  the 
slope  dW/dX  at  X]  (=  DWI). 
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SOLUTION  PROCEDURE 

First  we  guess  the 
Equations  [85]  and  [87] 
tlai  conditions  at  X  - 
and  Xm  are  determined 
tlons  in  Eqs.  [86]  and 
one  boundary  equation  1 
face.  With  the  Initial 
Eq .  [87]  at  X  *  X^,  we 
The  guessing  values  of 
until  Eq.  [91 ]  and  Eq. 


RESULT 


,alues  of  Xi  and  OKI. 
are  solved  with  the  1 n  1  - 
X] ,  and  the  values  of  X2 
with  extra  boundary  condl- 
[88].  In  principle,  only 
s  needed  for  each  lnter- 
condltlons  obtained  from 
solve  Eqs.  [89]  and  [90]. 
X]  and  DM  are  adjusted 
[78]  are  satisfied. 


The  typical  cavity  shapes  and  the  centerline 
pressure  distributions  are  shown  In  Figs.  16  and 
17.  Compared  with  the  numerical  schemes,  with  the 
current  model  we  only  need  to  guess  two  values 
rather  than  the  whole  cavity  shape,  and  the 
partial  differentall  equation  reduces  to  odlnary 
differential  equations.  A  large  amount  of  numeri¬ 
cal  work  is  saved  and  not  too  much  accuracy  is 
lost.  We  are  currently  obninlng  the  numerical  re¬ 
sults  by  finite-element  methods.  Compared  with 
the  short  bearing  model  (58),  the  present  model  ob¬ 
tains  better  results  when  correlated  with  the 
finite-element-method  model. 
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THE  INFLUENCE  OF  CAVITATION  ON  JOURNAL-BEARING  ROTORDYNAMIC  CHARACTERISTICS 
M.L.  ADAMS” 


Even  under  the  smoothest  running  conditions, 
the  journal  always  exhibits  some  orbital  vibration 
about  Its  equilibrium  position  with  respect  to  the 
bearing.  Provided  this  vibration  is  sufficiently 
small,  one  can  justify  characterization  cf  the 
resulting  fluid-film  Irteractive  dynamic  force 
between  journal  and  bearing  by  a  linear  mathemati¬ 
cal  model.  The  essential  Ingredients  of  this  model 
are  the  linearized  stiffness  and  damping  coeffici¬ 
ents  for  the  journal  bearing,  l.e.,  a  2  by  2  stiff¬ 
ness  matrix  and  a  2  by  2  damping  matrix.  These 
rotordynamic  coefficients  for  the  Interactive  dy¬ 
namic  force  between  the  journal  and  bearing  are  de¬ 
lineated  In  the  following  equation: 
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The  radial  dynamic  force  f  exerted  on  the  jour¬ 
nal  Is  the  vector  sum  of  the  applied  static  force 
W  and  the  total  fluid-film  force  F,  as  Illus¬ 
trated  In  Fig.  18. 


In  the  simplest  of  cases,  the  flow  field  be¬ 
tween  the  rotating  and  nonrotating  members  has  no 
memory.  This  would  be  the  situation  for  the  hydro- 
dynamic  pressure  distribution  produced  In  a  purely 
viscous,  fully  developed,  Incompressible  fluid- 
film  journal  bearing  in  the  absence  of  cavitation 
(film  •  iture)  as  characterized  by  the  following 
cl  ass  1 v  lubrication  equation  of  Reynolds  com¬ 
bined  with  pressure  boundary  conditions: 


<L/V  3fi\ 

3n\^i  dr\J 
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where  n  Is  the  unwrapped  circumferential  coordi¬ 
nate  and  5  is  the  axial  coordinate  map  of  the  cy¬ 
lindrical  bearing  surface. 

The  lubricant  fluid  film  pressure  distribu¬ 
tion  p(n,(3,t)  embodies  an  interactive  radial 
force  vector  which  acts  to  separate  the  journal 
and  the  bearing.  This  pressure  distribution  is 
functionally  related  to  the  radial-orbit  position 
and  velocity  vectors  of  the  journal  center  rela¬ 
tive  to  the  bearing  center;  that  Is,  the  film 
thickness  distribution  h  separating  the  journal 
and  the  bearing  is  a  function  of  the  relative  posi¬ 
tion  of  the  journal  with  respect  to  the  bearing  as 
follows: 


h  =  C  -  x  cos(r)/R>  -  y  s1n(n/R>  [94] 


1  department  of  Mechanical  and  Aerospace 
Engineering,  Case  Western  Reserve  University, 
Cleveland,  OH. 


where  x  *  xs  +  x  and  y  -  ys  +  y  (see  Fig.  18). 
Equation  [94]  Is  applicable  to  the  typically 
assumed  condition  that  bearing  clearance  Is  much 
smaller  than  bearing  radius  (C  «  R)  and  that  the 
bearing  and  the  journal  axial  centerlines  are  par¬ 
allel;  l.e.,  no  misalignment,  as  given  by  3h/3fj 
=■  0.  Consistent  with  Eq.  [94],  the  following 
derivatives  of  h,  which  appear  in  Eq.  [93],  are 


(a) 


W:  APPLIED  STATIC  LOAD 
F:  TOTAL  FLUID  FORCE 
f  :  DYNAMIC  FLUID  FORCE 


- »~y 

f  =  W  +  F 


(b) 


Fig.  18.  -  (a)  Radial  motion  coordinates  of  jour¬ 
nal  relative  to  bearing,  referenced  to  static 
equilibrium  state  (W  =  static  load);  (b)  Harmon¬ 
ic  orbital  motion  of  journal  relative  to  bearing; 
x  =  x0  sin(flt  +  0X);  y  =  y0  slnCQt  +  dy) . 
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written  as 

dh/dn  =  C(x  s1n(n/R>  -  y  cos(n/R)]/R  C95 j 

dh/dt  =  -x  cos(n/R)  -  y  sin(n/R)  C963 

Clearly,  In  this  case,  the  Instantaneous  pressure 
distribution  and,  thus,  the  Instantaneous  journal¬ 
bearing  separating  force  are  determined  by  the 
Instantaneous  position  and  velocity  vectors  (x,  y) 
and  (x,  y>,  respectively,  and  thus,  has  no  memory 
feature.  With  no  memory,  the  eight  stiffness  and 
damping  coefficients  are  functions  only  of  the 
static  equilibrium  flow  field,  e.g.,  for  purely 
laminar  flow,  a  function  of  bearing  Sommerfeld 
number  only. 

When  cavitation  Is  admitted  Into  the  problem, 
the  longstanding  approach  Is  to  require  p  =  pvap 
and  9p  •  h  =  0  (to  conserve  mass  flow)  at  the 
Interface  boundary  separating  the  full-film  region, 
In  which  Eq.  [93]  Is  solved,  and  the  cavltated  re¬ 
gion,  In  which  all  pressures  are  set  equal  to  the 
vapor  pressure,  pvap.  With  this  approach  to  cavi- 
tatlng  bearings,  the  solution  again  leads  to  rotor- 
dynamic  coefficients  that  are  functions  of  the 
static  equilibrium  flow  field  only.  While  this  ap¬ 
proach  to  cavltatlng  bearings  Is  probably  an  ac¬ 
ceptable  characterization  for  static  load  capacity 
predictions,  It  has,  unfortunately,  been  extended, 
with  little  questioning,  to  the  rotor-dynamic  coef¬ 


ficients.  If  closely  scrutinized,  the  approach 
for  cavltatlng  bearings  is  found  to  satisfy  mass 
conservation  only  within  the  full-film  region  and 
on  Its  boundaries  but  not  within  the  cavltated  or 
ruptured  region.  In  essence,  this  approach  func¬ 
tions  within  a  rotordynamlc  context  In  the  follow¬ 
ing  loosely  worded  manner:  whenever  the  right-hand 
side  of  the  Reynolds  equation  gives  the  command, 
the  left-hand  side  Instantly  becomes  activated  to 
allow  the  generation  of  local  film  pressures, 
whether  or  not  the  fluid  has  had  time  to  refill 
the  local  gap. 

In  fact,  the  actual  system  Is  quite  transient 
or  memory-like,  as  shown  by  the  computations  of 
Elrod  and  Adams  (14),  who  employed  a  time-transient 
solution  scheme  which  properly  accounted  for  mass 
conservation  throughout  tne  bearing  clearance  gap. 
Thus,  even  In  the  absence  of  fluid  inertia,  a  dy¬ 
namically  loaded  cavitating  journal  bearing  can  ex¬ 
hibit  Inertia-like  effects.  And  these  effects  can 
be  expected  to  become  progressively  more  important 
as  vibration  frequency  is  increased.  This  is 
quite  similar  to  compressibility  effects. 

A  more  complete  discussion  of  this  and  related 
topics  is  given  In  (59).  For  the  purpose  of  this 
discussion,  it  Is  simply  concluded  that  the  memory 
or  frequency  effects  of  cavitation  on  bearing 
rotordynamlc  characteristics  is  Important. 
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CAVITATION  IN  JOURNAL  BEARINGS  AND  SQUEEZE-FILM  DAMPERS:  ANALYTICAL 


B.  JACOBSON14 

Cavitation  boundary  conditions  for  stationary 
and  vibrating  journal  bearings  have  been  discussed 
In  the  literature  for  more  than  30  years.  The 
Swlft-Stleber  condition  for  the  start  of  the  cavi¬ 
tation  region  has  been  used  both  for  stationary 
bearings  and  whirling  bearings  where  the  minimum 
film  thickness  Is  constant  and  the  speed  of  the 
cavitation  zone  Is  less  than  the  mean  value  of  the 
surface  speeds.  A  more  accurate  boundary  condi¬ 
tion  for  the  start  of  the  cavitation  zone  was 
given  by  Floberg  where  the  number  of  streamers  per 
unit  width  In  the  cavitation  zone  determined  the 
amount  of  subcavity  pressure  In  front  of  the  cavi¬ 
tation  zone.  The  most  frequently  used  outlet 
boundary  condition  for  the  cavitation  zone  Is  the 
JFO  < Jal-  3bs son ,  Floberg,  Olsson)  condition,  which 
says  that  the  oil  flow  going  Into  the  cavitation 
region  Is  transported  straight  through  It  without 
side  flow  and  that  there  Is  continuity  of  mass 
flow  at  the  end  of  the  cavitation  region. 

THE  DIFFERENCE  BETWEEN  VAPOR  AND  GAS  CAVITATION 

Whereas  the  pressure  In  a  vapor  bubble  Is 
almost  totally  given  by  the  temperature,  the 
pressure  In  a  gas  bubble  Increases  rapidly  accord¬ 
ing  to  the  gas  law  as  the  volume  of  the  bubble 
decreases.  The  transport  of  gas  Into  and  away 
from  a  bubble  Is  a  slow  process,  as  can  be  clearly 
seen  from  the  cavitation  films  where  It  often 
takes  hundreds  of  vibrations  before  the  gas  cavita¬ 
tion  pattern  stabilizes.  If  the  gas  transport 
Into  and  out  of  the  bubble  is  neglected,  the  bub¬ 
ble  pressure  can  easily  be  calculated  using  a  poly- 
trop'c  compression,  l.e.,  midway  between  an  Iso¬ 
thermal  and  an  adiabatic  compression.  This  gives 
a  radius  ratio  of  0.2  for  a  flat  bubble  when  Its 
pressure  Is  Increased  100  times.  Therefore,  If 
the  bubble  contains  air,  It  will  be  visible  In  the 
high-pressure  zone  of  the  bearing.  Vapor  bubbles, 
on  the  other  hand,  collapse  within  milliseconds  of 
formation. 

ANALYSIS  OF  VAPOR  CAVITATION 


A  simple  model  of  the  vapor  cavitation  behav¬ 
ior  In  dynamically  loaded  journal  bearings  Is  to 
assume  the  pressure  to  be  constant  during  the  col¬ 
lapse  of  the  bubble.  This  means  that  the  vapor 
condenses  at  the  same  rate  as  the  bubble  wall 
moves.  For  a  spherical  bubble,  this  gives  the 
collapse  velocity  as 

MTV  , 

V  =  C97] 

V 

and  for  a  flat  bubble  with  the  radius  r  In  an 
oil  film  of  the  thickness  h, 


1 ’Chalmers  University  of  Technology,  Gote- 
borg,  Sweden,  and  SKF  Engineering  and  Research  Cen¬ 
tre,  Nleuwegeln,  The  Netherlands. 
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[98] 


where  hc  Is  the  film  coefficient  of  heat  trans¬ 
fer  for  condensing  oil  vapor  at  low  flow  velocity, 
/v  Is  the  heat  of  evaporation  for  oil,  and  ATvj 
Is  the  temperature  drop  from  vapor  to  liquid  oil! 

For  a  10-mm-d lame  ter,  flat  bubble  In  a  20-pm 
oil  film  and  with  p  -  6xl0~6  kg/m3,  he  »  600 
W/m2*°C,  Zv  «  350  kJ/kg,  and  the  collapse  velocity 
V  -  5  m/s,  the  temnorature  drop  needed  from  the  va¬ 
por  to  the  oil  «  j  »  7.0x1 0-3  C.  This  means 
that  the  vapor  t „s  collapse  almost  without  any 
resistance. 


CAVITATION  FORMATION 


By  changing  the  coordinate  system  for  the 
vibrating  bearing  to  a  system  rotating  with  the 
same  speed  as  the  minimum  oil  film  thickness,  the 
right-hand  side  of  the  Reynolds  equation  can  be 
simplified  for  times  close  to  the  Instant  the  mini¬ 
mum  oil  film  thickness  has  Its  lowest  value.  The 
right-hand  side  of  the  Reynolds  equation  for  the 
geometry  of  the  vibrating  bearing  Is  then  simpli¬ 
fied  to 


sin  0 


Vcs  +  cd) 


“  cdud 


[99] 


For  close  to  zero  minimum  oil  film  thickness,  es  + 
cj  *  1 ,  and  the  criterion  for  the  formation  of 
vapor  cavitation  Is  simply 


cd 


[100] 


Just  as  was  mentioned  In  the  experimental  results, 
a  high  rotating  speed  ws  of  the  journal  stops 
the  vapor  cavitation  (see  also  (24)).  This  Is 
easily  seen  as  ed  must  always  be  less  than  1. 


CAVITATION  COLLAPSE 


Once  vapor  cavitation  has  been  formed  and 
some  time  has  passed,  the  minimum  oil  film  thick¬ 
ness  Increases,  and  the  atmospheric  pressure  at 
the  edges  of  the  bearing  becomes  Important  for  con¬ 
sideration  of  the  collapse  of  the  vapor  cavitation 
bubbles.  The  side  flow  Induced  by  the  atmospheric 
pressure  at  the  bearing  edges  increases  with  the 
third  power  of  the  oil  film  thickness.  The  high¬ 
speed  films  showed  that  the  bubble  collapse  Is 
mainly  governed  by  this  side  flow.  A  simple  short 
bearing  type  of  analysis  can  give  the  criterion 
for  when  the  bubble  reaches  its  maximum  width  and 
starts  to  collapse.  This  happens  when 
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(d)  BUBBLE  EXPANDS  (t  =  R2.51  MS;  <e)  BEGIN  SEPARATION  (t  =  33.89  MS;  (f)  OCCURRENCE  OF  PEAK  PRESSURE 

£  =  0.83).  £  =  0.90).  t  =  29.31  MS;  e  =  G.88). 


Fig.  19.  -  Pressure  distribution  and  bearing  configuration  for  a  full  period 
of  shaft  whirl  (0.1  s  e  s  0.9).  Figures  (a)  to  (f)  viewed  clockwise  are  con¬ 
secutive  in  time. 


(Ua  +  V  3h  3h  h3  <po  "  pc}  [101] 

2  ax  +  9t  6n  ,  ) 2 

m  “  zc 

If  the  cavitation  collapse  takes  place  in  the  same 
location  many  times,  It  will  result  in  cavitation 
erosion  damage  on  the  journal  and  the  bearing. 


NUMERICAL-THEORETICAL  CALCULATIONS  OF 
VAPOR  CAVITATION 

The  appearance  and  motion  of  the  cavitation 
zone  and  its  boundaries  can  be  studied  using  numer¬ 
ical  methods.  This  was  done  using  a  computer  pro¬ 
gram  (7)  in  which  Elrod's  algorithm  for  the  conti¬ 
nuity  of  mass  flow  into,  through,  and  out  of  the 
vapor  cavitation  zone  was  used.  This  program  makes 


it  possible  to  predict  the  dynamic  behavior  of  the 
vapor  cavitation  zone.  The  Cray  computer  at  NASA 
Lewis  Research  Center  was  used  to  do  the  computa¬ 
tions.  The  boundary  condition  for  the  vapor  cavi¬ 
tation  zone  was  the  pressure  equal  to  zero,  which 
is  consistent  with  the  earlier  experimental  re¬ 
sults.  This  was  used  together  with  Elrod's  algo¬ 
rithm  to  make  sure  that  the  continuity  of  mass 
flow  was  maintained,  both  in  the  full  oil  film  and 
In  the  cavitation  zone,  during  the  whole  vibration 
cycle.  Some  numerically  calculated  pressure  dis¬ 
tributions  and  cavitation  regions  (23)  are  shown 
In  Fig.  19.  A  comparison  with  the  experimental  re¬ 
sults  shows  similar  behavior  for  the  experimental 
cavitation  zone  and  the  theoretically  calculated 
zone,  even  if  the  size  of  the  experimental  cavita¬ 
tion  zone  is  slightly  smaller.  This  can  be  partly 
explained  by  elastic  distortion  of  the  PMMA  tube. 
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CONCLUSIONS 

Vapor  cavitation  occurs  when  the  tensile 
stress  applied  to  the  oil  exceeds  the  tensile 
strength  of  the  oil  or  the  binding  of  the  oil  to 
the  surface.  The  physical  situation  necessary  for 
vapor  cavitation  to  exist  is  a  squeezing  and  slid¬ 
ing  motion  In  a  bearing.  The  Implication  of  the 
results  as  they  apply  to  squeeze  film  dampers  is 
that  if  the  sliding  velocity  Is  increased,  the 
occurrence  of  vapor  cavitation  is  decreased  (see 
also  (24)). 


Vapor  cavitation  bubbles  collapse  at  pressures 
lower  than  the  atmospheric  pressure,  and  they  can¬ 
not  be  transported  through  a  high  pressure  zone  of 
the  bearing,  unlike  gas  cavitation  bubbles. 

The  numerical  calculations  (23)  of  the  cavita¬ 
tion  bubble  dynamics  show  that  an  increasing  vibra¬ 
tion  amplitude  gives: 

•  Longer  vapor  bubble  life 

•  Larger  maximum  vapor  bubble  size 

•  Earlier  start  and  later  collapse  of  the  va¬ 
por  bubble. 
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The  following  discussion  summarizes  the 
major  issues  brought  out  by  the  contributors  as 
these  issues  relate  to  present  technological 
requirements.  At  the  risk  of  appearing  myopic, 
the  editors  have  taken  the  liberty  to  occasion¬ 
ally  emphasize  the  importance  of  the  work,  what 
has  been  resolved,  and  what  needs  to  be  done. 

Many  of  the  articles  focused  on  cavitation 
as  it  pertained  to  a  broad  list  of  applications 
in  crankshaft  bearings,  rotordynamics,  slider 
bearings,  and  seals.  The  key  issue  is  the 
dynamic  loading.  As  pointed  out  in  the  intro¬ 
duction  to  Part  II,  dynamic  loading  causes 
changes  in  the  local  film  thickness,  which  leads 
to  nonstationary  cavitation.  This  certainly 
plays  an  important  role  in  determining  rotordy- 
namic  characteristics.  Adams'  discussion  ampli¬ 
fies  this  topic,  especially  concerning  the 
effect  of  frequency.  Adams  draws  the  analogy 
that  nonstationary  cavitation  can  exhibit 
"inertia-like  effects"  and  that  these  effects 
can  be  expected  to  become  progressively  more 
important  as  vibration  frequency  is  increased. 

Of  course,  this  is  only  one  example  of  many 
where  nonstationary  cavitation  is  likely  to  have 
a  significant  effect  on  the  rotordynamic  charac¬ 
teristics.  Jacobson  discusses  the  effect  of 
deformation,  vibrational  amplitude,  surface  ten¬ 
sion,  and  sliding  velocity  on  vapor  bubble 
inception,  size,  life,  and  collapse.  Both 
inception  and  collapse  play  an  important  role  in 
the  bubble  dynamics  which,  in  turn,  affects 
rotordynamic  behavior. 

As  pointed  out  by  Launay,  dynamic  loading 
is  the  primary  concern  in  the  design  of  journal 
bearings  of  automotive  engines.  Here,  operation 
with  an  incomplete  film  is  typical  in  such  bear¬ 
ings,  and  one  strives  to  maintain  an  adequate 
film  while  keeping  the  friction  as  low  as  pos¬ 
sible.  In  the  case  of  diesel  engines,  the  com¬ 
bustion  conditions  can  become  more  severe  (i.e., 
higher  rates  of  change  of  cylinder  pressure), 
resulting  in  cavitation  erosion  damage  to  the 
bearings.  Martin  discusses  the  work  on  this 
subject  by  Garner,  James,  and  Warriner  (52).  He 
Includes  examples  of  typical  cavitation  erosion 
damage,  together  with  a  plausible  mechanism  of 
damage.  Commonly  applied  palliatives  and  their 
effectiveness,  based  on  engine  experience,  are 
also  discussed.  It  was  concluded  that  cavita¬ 
tion  erosion  damage  resulted  from  three  main 
causes.  Two  were  related  to  oil-column  inertia 
effects  in  the  connecting  rod  drillings,  crank¬ 
shaft  drillings,  and  bearing  oil  feed  features. 
These  were  obviated  by  simple  modifications  to 
the  detail  design.  The  third,  related  to  jour¬ 
nal  movements  within  the  beaming  clearance 
space,  posed  a  more  difficult  problem.  Martin 
indicated  that  a  more  detailed  theory  than  that 
used  in  their  analysis  was  needed  to  reliably 
predict  damage  for  the  design  stage. 

Of  course,  cavitation  effects  are  not  lim¬ 
ited  to  journal  bearings  and  squeeze-film  damp¬ 
ers.  In  a  seal  or  slider  application,  Lebeck 


describes  a  mechanism  in  which  cavitation  is 
used  to  advantage.  The  asperities  in  parallel 
sliders  were  hypothesized  to  produced  load  sup¬ 
port  by  virtue  of  the  fact  that  cavitation  deve¬ 
loped  at  the  downstream  side  of  the  asperity. 

In  the  presence  of  cavitation  there  is  a  net 
load  support  for  individual  asperities  and, 
hence,  for  the  parallel  slider. 

The  above  articles  dealt  mainly  with  cavi¬ 
tation  effects  in  various  applications,  especially 
where  it  concerned  dynamic  operating  conditions. 

In  order  to  gain  Insight  and  provide  better 
design  through  numerical  modeling  for  these 
applications,  we  must  identify  the  modes  of  cav¬ 
itation  associated  with  different  dynamical  con¬ 
ditions.  For  example,  in  the  introduction  to 
Part  I  we  distinguished  between  gaseous  and 
vaporous  cavitation,  but  in  actual  practice 
both  may  occur  together  unless  definitive  steps 
are  taken  to  preclude  one  or  the  other  (i.e., 
degassing  the  fluid  so  that  only  vapor  cavita¬ 
tion  can  occur  or  pressurizing  the  fluid  so  that 
vapor  pressures  are  not  possible).  Jacobson 
elaborates  on  some  of  the  noted  differences 
between  vapor  and  gaseous  cavitation  in  Part  II. 
Parkins'  article  on  the  "Categorization  of  Cavi¬ 
tation  Types  and  Onset  Conditions"  is  a  begin¬ 
ning  in  understanding  the  conditions  under  which 
one  or  the  other  might  dominate.  The  categori¬ 
zations  in  this  article  were  based  on  sustained 
oscillatory  squeeze  motion  -  simplified  by 
removing  the  translational  component.  Three 
cavitation  regimes  -  including  subregimes  -  were 
Identified  and  characterized  according  to  the 
method  of  feeding  the  bubble.  It  was  noted  that 
cavitation  bubbles  fed  from  withi::  the  film 
could  be  either  gaseous  or  vapor.  Cavitation 
bubbles  fed  from  outside  the  film  obviously  were 
gaseous  (e.g.,  refer  to  the  proposition  of  "side 
entrainment  of  gas,"  by  Heshmat)  and  took  sev¬ 
eral  forms.  Three  of  these  forms  were  leaf  type 
patterns,  and  the  fourth  appeared  as  smoothly 
curved  bubbles  entering  from  all  sides  of  the 
square  plate.  Certain  regimes  were  observed  to 
coexist.  Parkins  Indicates  that  "a  limited 
knowledge  is  available  about  the  conditions  which 
determine  the  onset  of  a  particular  regime  or 
subregime  and  its  transformation  into  another 
subregime.  The  large  number  of  possible  forms 
of  behavior  make  a  confused  picture  when 
observed  over  a  small  range  of  variables,  which 
renders  it  difficult  to  formulate  classifica¬ 
tions."  In  the  editors'  opinions,  this  work 
deserves  an  expanded  effort  in  which  classifica¬ 
tion  according  to  methods  of  feeding  the  bubble 
should  be  looked  at  with  the  translational  com¬ 
ponent  as  well.  It  is  only  through  the  classi¬ 
fication  procedure  that  we  can  hope  to  bring 
order  to  this  "confused  picture,"  thereby  pro¬ 
viding  a  basis  upon  which  analysts  can  formu¬ 
late  appropriate  models. 

Me  have  already  noted  the  importance  of 
bubble  dynamics  in  rotordynamic  behavior.  Bub¬ 
ble  dynamics  is  a  complicated  subject  in  itself. 
The  dynamics  of  the  bubble  can  be  affected  by 
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the  gas  and  vapor  content  (or  fractional  film 
content  within  the  bubble)  as  well  as  by  trans¬ 
cavity  flow  rate  or  fluid  propagation  speed 
within  the  bubble.  The  oropagation  speed  of  the 
fluid  through  the  bubble  depends  on  the  morpho¬ 
logical  condition  of  the  bubble.  Pan  and  Frene 
discuss  this  issue  with  regard  to  numerical 
modeling  and  implementing  the  transport  process 
across  the  interface.  They  include  the  possi¬ 
bilities  that  a  portion  of  the  void-dominated 
region  is  locally  filled  with  liquid,  oridging 
both  surfaces,  or  that  a  thin  liquid  layer  may 
be  attached  to  either  but  not  both  surfaces 
(also  see  discussions  by  Heshmat  by  Barrett  and 
Branagan  in  Part  I).  For  the  first  possibility, 
the  propagation  speed  of  the  fluid  is  assumed 
to  be  the  average  of  the  two  surface  speeds. 

But  for  the  latter  possibility,  the  fluid  takes 
on  the  speed  of  the  surface  to  which  it  becomes 
attached.  In  the  case  of  the  journal  bearing, 
one  might  argue  that  fluid  transport  through  the 
cavity  influences  the  position,  size,  and  life 
of  the  bubble.  Consequently,  it  can  have  an 
effect  on  the  load  capacity  (7)  and,  possibly, 
the  stiffness  and  damping  characteristics  of  the 
bearing.  Knowledge  is  scant  on  the  conditions 
that  determine  whether  attachment  is  to  one  or 
both  surfaces,  although  it  is  generally  accepted 
that  load  is  a  determining  factor.  More  work 
is  needed  to  clarify  this  issue  and  the  effect 
it  has  on  bubble  dynamics  and,  consequently,  on 
rotordynamic  performance. 


The  temperature  effects  on  cavitation  are 
an  Important  consideration  both  from  the  stand¬ 
point  of  bubble  dynamics  and  heat-transfer 
characteristics.  Both  Heshmat  and  Barrett  and 
Branagan  discussed  temperature  effects  as  it 
pertained  to  "hot-oil  carry  over"  in  padded 
bearings.  Not  surprisingly,  it  was  noted  by 
both  that  heat  generation  within  the  cavitation 
region  was  practically  nonexistant.  Barrett  and 
Branagan  compared  theory  with  experiment.  They 
found  that  treating  the  void  region  as  a  two- 
phase  lubricant-gas  in  which  the  fluid  attaches 
Itself  to  the  shaft  was  the  most  suitable  model 
for  the  experiment  considered.  Further  studies 
involving  a  wider  range  of  operating  conditions 
would  be  useful  information  to  the  analyst  in 
modeling  thermal  effects  on  cavitation.  Cor¬ 
rectly  predicting  the  apportionment  of  heat 
transfer  to  the  shaft  and  housing  is  partly 
dependent  on  choosing  the  appropriate  cavita¬ 
tion  model.  One  other  way  in  which  thermal 
effects  can  affect  cavitation  is  through  the 
Indirect  changes  in  surface  energy.  In  this 
context  we  refer  to  temperature  effects  on  the 
wettability  of  the  lubricant  to  the  surface  and 
the  consequences  to  fluid  transport  through  the 
bubble  as  well  as  heat  transfer  to  the  bearing 
and  shaft  material.  Albeit,  Szeri,  Rajagapol, 
and  Zhang  (Part  II)  noted  that  journal  bearings 
lubricated  with  non-Newtonian  lubricants  could 
be  less  susceptible  to  thermal  effects  when  com¬ 
pared  with  bearings  lubricated  with  Newtonian 
fluids. 


Non-Newtonian  effects  on  cavitation  were 
further  discussed  by  Szeri  and  shown  to  be 
important.  That  is,  a  shear  thickening  fluid 
cavitated  within  a  range  of  material  parameters, 
whereas  the  shear  thinning  fluid  did  not  cavi¬ 
tate  at  all  for  the  same  parameter  range. 
Furthermore,  the  change  in  material  parameter 
for  the  shear  thickening  fluid  affected  the 
position  of  th  rupture  boundary  to  some  extent. 
Since  many  of  the  lubricants  used  these  days  are 
non-Newtonian,  it  behooves  us  to  learn  much 
more  about  this  phenomena. 

Using  high-speed  photography,  Jacobson  and 
Hamrock  (Part  I)  noted  a  clear  influence  of  the 
surface  energy  of  the  journal  surfaces  on  cavi¬ 
tation  incipience.  For  steel  journals,  cavita¬ 
tion  bubbles  start  in  different  parts  of  the 
low-pressure  region,  depending  on  where  the 
largest  bubble  or  dirt  particle  is  located.  For 
PTFE  journals,  the  cavitation  seems  to  start 
from  the  same  position  of  the  journal  for  each 
vibration,  which  Indicates  a  rupture  of  the  oil- 
PTFE  interface.  Presently,  analytical  models  do 
not  account  for  this  phenomena. 

This  calls  attention  to  the  subject  of 
numerical  modeling.  Most  of  the  contributions 
to  this  subject  have  dealt  with  flow  preserving 
schemes,  numerical  expedience,  and  some  recent 
advances  using  finite-element  schemes.  In  flow 
preserving  schemes,  the  attention  is  naturally 
drawn  to  the  more  complex  cavitation  boundaries. 
With  the  advent  of  the  Elrod  algorithm,  much  of 
the  complexity  has  been  reduced  while  retaining 
the  flow-preserving  features  of  the  JFO  theory. 
Keith  and  VI jayaraghavan  (Part  II)  present  their 
version  of  the  Elrod  algorithm  in  the  context  of 
gas  dynamics.  Borrowing  concepts  used  in  model¬ 
ing  transonic  flow  through  a  duct,  an  artificial 
viscosity  term  is  introduced  that  implements  the 
"type  differencing"  employed.  The  two  schemes 
are  equivalent,  but  the  artificial  viscosity 
plays  an  opposite  role  to  Elrod's  switch  func¬ 
tion.  Their  formulation  was  instructive  and  not 
unlike  the  modeling  viewpoint  of  Pan  and  Frene 
(see  discussion  on  Rankine-Hugoniot  conditions 
in  Part  II),  although  Pan  and  Frene's  treatment 
is  more  general  in  that  it  accomodates  any  one 
of  three  types  of  fluid  transport  within  the 
cavitation  boundary.  Frene,  Launay,  and  Pan 
illustrate  the  use  of  the  equivalent  Rankine- 
Hugoniot  conditions  for  a  centrally  fed  journal 
bearing  using  short  bearing  theory.  The  analogy 
drawn  between  the  shock  wave  and  the  cavitation 
interface  is  constructive  since  it  forms  a  basis 
of  communication  and  progress  in  either  field 
can  be  shared. 

As  our  models  become  more  complex,  the 
search  for  simpler  more  efficient  programing 
techniques  becomes  essential.  Tichy  and  Ku 
(Part  II)  describe  a  solution  procedure  for  a 
squeeze-film  damper  bearing  with  two  circumfer¬ 
ential  oil  supply  grooves  at  the  bearing  ends 
undergoing  circular-centered  whirl.  In  this 


64 


procedure  a  large  amount  of  numerical  worK  is 
saved  and  not  too  much  accuracy  is  lost.  Numer¬ 
ical  results  for  centerline  pressure  and  cavity 
shape  using  finite-element  methods  compared 
favorably. 

This  brings  us  to  our  last  topic  of  discus¬ 
sion,  employing  finite-element  methods  in  lubri¬ 
cation  problems  in  which  cavitation  occurs. 
Because  of  the  ease  in  programing  and  adaptabil¬ 
ity  to  the  high-speed  super  computers,  finite- 
difference  schemes  are  most  often  used.  How¬ 
ever,  from  the  standpoint  of  handling  irregular 
boundaries,  the  finite-element  method  is  well 
suited  to  dealing  with  cavitation.  In  Part  II 
Booker  discusses  several  mathematical  models 
commonly  used  to  incorporate  cavitation  in  the 
finite-element  analysis  of  lubrication.  The 
models  are  classic  in  that  they  are  not  flow¬ 
preserving.  Until  recently,  this  has  been  an 
obstacle  to  achieving  a  physically  realistic 
solution  to  transient  problems.  Bayada  (Part 
II)  discusses  a  two-dimensional  analysis  using 
FEM  in  which  mass  conservation  is  treated  in  a 
way  that  is  consistent  with  that  proposed  by 
Elrod  (14)  or  the  equivalent  Rankine-Hugoniot 
condition  at  the  interface.  The  procedure  was 
applied  to  a  seal  (46)  using  rectangular  finite 


elements.  The  discrepancy  between  input  mass 
flow  and  output  mass  flow  never  exceeded 
4  percent  in  each  period.  As  of  this  writing, 
Booker  (59)  has  achieved  a  flow-preserving  FEM 
scheme  applied  to  journal  bearings  in  three 
dimensions. 

To  sum  up,  the  experimenter  must  continue 
to  help  perceive  the  order  in  a  very  complex 
phenomena.  This  can  be  done  by  extensive  clas¬ 
sification  and  clarification  of  unresolved 
issues  using  carefully  planned  experiments.  The 
numerical  analyst  find  themselves  on  a  tortuous 
journey  to  realism  in  which  all  of  the  effects 
of  temperature,  misalignment,  deflections,  local 
deformations,  surface  energy,  and  lubricant  rhe¬ 
ology,  to  name  only  a  few,  must  eventually  be 
included  in  the  computational  models.  To  make 
this  possible,  the  analysts  must  bring  to  bear 
the  latest  numerical  advances,  keeping  abreast 
of  new  developments  in  related  fields.  Algo¬ 
rithms  that  incorporate  all  of  the  latest  tech¬ 
niques  in  finite  differencing  as  well  as  finite 
elements  to  shorten  computational  times  are 
needed.  In  addition,  these  same  algorithms  must 
lend  themselves  to  employing  the  latest  computer 
architecture,  i . e . ,  pipeline  and  parallel 
processing. 
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